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ABSTRACT 

The research objective was to measure near-angle 
scattering (i.e., scattering less than one degree from the 
specular beam) from a binary optic and determine if a scalar 
scattering theory could adequately predict scattering 
performance. Near-angle scattering was measured from a binary 
optic beam splitter. A scalar scattering theory was developed 
and modeled using FORTRAN on a personal computer; results from 
the computer simulation are compared to the actual 
measurements. 

The scalar scattering theory modifies Fraunhofer 
diffraction by including two types of surface topography which 
contribute to the scattering, specifically: 1) surface 

roughness due to micro-irregularities which are considered to 
be randomly distributed and 2) large scale surface features 
(i.e., the binary step pattern) which are deterministic. The 
random surface roughness, the autocorrelation length of the 
roughness, and the height of the binary optic's pattern were 
determined using a Talystep surface prof ilometer . The scalar 
theory appears to give good results when compared to the 
measurements. However, higher values were required for the 
surface roughness and autocorrelation lengths than the 
Talystep indicated, which may be due to the Talystep' s short 
scan length. 



iv 



X - 



Military and space applications for binary optics are also 
discussed . 
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I. 



INTRODUCTION 



A. BINARY OPTICS 

Binary optics is a revolutionary optical technology. 
Quality conventional optics are dependent on the skill of the 
lensmaker and the geometric form of the optic's surface. 
However, binary optics are designed with the aid of a computer 
and very large scale integrated (VLSI) lithographic 
technology. Binary optics get their name from the binary 
profile etched onto the optic's surface (Figure 1.1). 




Binary optics is also called diffractive optics. However, 
diffractive optics commonly refers to a multi-level or 
"staircase" pattern on the surface of the optic. The 

staircase pattern (Figure 1.2) is achieved by repeated 
application of binary masks, which is done to improve the 
efficiency. 
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Figure 1.2 Profile of Diffractive (Binary) Optic 



1 . Definition of Binary Optics 

Wilfrid Veldkamp, the developer of binary-optics 
technology, defines binary optics as: 

Binary optics is a diffractive optics technology that 
uses computer generated designs of microscopic relief 
patterns and VLSI ion-etching technology. It creates 
novel diffractive optical devices and provides freedom and 
new materials choices for conventional refractive 
elements. It is a classical example of a technology 
transfer — from electronics to optics (and back) . Through 
advances in VLSI technology patterns with characteristic 
dimensions as small as the wavelength of light can now be 
mass produced. The name binary optics, was coined to 
reflect a surface relief optics technology based on (1) 
a binary electronic mask technology (chrome on glass) , (2) 
a binary mask design code, 2 N phase steps for N masks, and 
(3) a two level or high-low, ion etch process with uniform 
surface treatment. [Ref. 1] 

2 . Advantages of Binary Optics 

Binary optics have several advantages over 
conventional optics, specifically: 

- Multi-level binary optics can be made highly efficient. 

- In a system design, use of binary elements may reduce the 
total number of elements, thereby saving weight. 
Additionally, certain applications may allow the binary 



2 



element to be embossed in plastic, which would provide a 
weight savings as well as the potential for cost savings. 

- Binary optics can perform functions that conventional 
elements cannot perform. 

- Binary optics allow flexibility in optical system design 
to generate compact optical shapes which will produce the 
desired optical wave front. 

- Binary optics can save cost by correcting aberrations in 
conventional lenses, allowing a lower quality conventional 
lens to be used in the system. Additionally, binary 
optics can be used with conventional lenses to form an 
equivalent asphere or achromatized element which will be 
at a lower cost than using conventional optics. 

B. APPLICATIONS FOR BINARY OPTICS 

Binary optics are very versatile; they can serve as thin 
lenses, prisms, gratings, holograms, or phase plates. 
Impressively, binary optics are able to perform functions 
impossible with conventional optics, specifically they can 
function as high speed rotary scanners, multiplexers, filters, 
laser beam shapers and combiners. Currently, lithographic 
fabrication resolutions primarily limit the use of binary 
optics to the infrared (IR) and visible wavelengths. Binary 
lenses can be made in mass quantity and thus have the 
potential to lower optical component costs. Additionally, the 
binary elements can replace heavier conventional elements and 
in some cases reduce the number of lenses required in an 
optical system, simplifying and compacting the design while 
reducing the weight. Binary optics technology has broad 
applicability to military and space systems. [Ref. 2] 
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Unfortunately, binary optics is not the panacea for all 
optical needs; binary optics will never totally replace 
conventional optics. Binary optics operate using diffraction 
and are limited to narrow bandwidth applications. However, 
when you combine binary optics with conventional optics many 
desirable results may be obtained (e.g., improved optical 
performance for lower cost) . 

C. SCATTERING ANALYSIS 

In the design and fabrication of optical systems, it is 
important to measure or know the level of scattering from 
optics. Experience indicates that even the best optical 
surfaces, fabricated using the best state-of-the-art (SOTA) 
polishing techniques, retain a microscopic scale surface 
roughness. Furthermore, even the best optical glasses have 
small local variations in their refractive indices. Both 
micro surface roughness and inhomogeneities within the glass 
produce refractive scatter on transmission. Additionally, 
micro surface roughness produces reflective scatter off the 
incident face of the optic, but this thesis addresses only the 
refractive scatter. 

1. Benefits of Scattering Analysis 

In a recent article concerning stray light analysis 
(scattering analysis falls under the definition of stray light 
analysis) , system level managers are strongly urged to "first 
perform the required analyses and then plan on measuring the 
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assembled system to the best of their capability." 
article also discusses benefits of scattering analysis which 
may be realized: [Ref. 3] 

- Enhance or change design. 

- Save misdirection of effort by preventing needless 
enhancements. 

- Support a major decision during the development of a 
sensor. 

- Perform tradeoff studies which can be used to help system 
level decisions. 

D. NEAR-ANGLE SCATTERING FROM A BINARY OPTIC 

1 . Near-Angle Scattering Measurements 

The research objective of this thesis was to measure 
near-angle scattering (i.e., scattering less than one degree 
from the specular beam) from a binary optic and determine if 
a scalar scattering theory could adequately predict scattering 
performance. Near-angle scattering is measured from a binary 
optic beam splitter manufactured by M.I.T. Lincoln Labs (see 
Chapter IV) . The measurements were taken using a HeNe laser 
which is part of a near-angle scattering instrument developed 
by Douglas W. Ricks at the Naval Weapons Center, China Lake, 
California. [Ref. 4] From the scattering measurements, a 
qualitative indication of fabrication errors and manufacturing 
limitations as applied to binary optics may be obtained. 

2 . Scalar Scattering Theory 

A scalar scattering theory was developed and modeled 
using FORTRAN on an IBM personal computer. Results from the 
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computer simulations were compared against the experimental 
scattering measurements (see Chapter V) . 

The scalar scattering theory modifies Fraunhofer 
diffraction by including two types of surface topography which 
contribute to the scattering, specifically: 

- Surface roughness from micro-irregularities are considered 
to be randomly distributed on the surface. 

- Large scale surface features (i.e., the binary step 
pattern) are deterministic. 

The theory requires knowledge of the binary optic's 
random surface roughness, the autocorrelation length of the 
roughness, and the height of the binary pattern, which were 
determined using a Talystep surface prof ilometer . In Chapter 
V, the model's surface roughness and surface roughness 
autocorrelation length were adjusted to obtain the best fit 
between the scalar theory and the measurements. The scalar 
theory appears to give good results when higher values are 
used for surface roughness and autocorrelation length than 
the Talystep measurements indicated. 

3 . Advantages of Using a Scalar Scattering Theory 

The primary advantage of using a scalar theory vice 
a full electromagnetic (EM) theory is in the relative 
simplicity of the theory and the reduction in computer time 
to perform the modeling, resulting in significant cost 
savings. A scalar theory is much simpler than a vector theory 
which must account for the phase and polarization of the 
light. In this thesis, a scalar scattering theory is 
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developed which is easily implemented on an IBM personal 
computer, thus providing an economical method of obtaining 
approximate performance based on the theory presented in 
Chapter V. The primary advantages of a scalar theory and the 
associated computer modeling, compared to the full EM theory, 
is twofold: 

- Full EM theory requires a highly knowledgeable 
experimenter, it is complicated and cumbersome. 

- Full EM theory is highly computer time intensive and 
therefore costly compared to a simple scalar model run on 
a personal computer (i.e., full EM theory may cost 
$100,000.00 using a full EM analysis program; alternately, 
the performance can be approximated on a personal computer 
for a few dollars using a scalar theory) . However, the 
scalar theory does not yield an exact solution but 
provides a good approximation. 
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II. FABRICATION OF BINARY OPTICS 



Binary optics is an emerging technology made possible by 
computer-aided design (CAD) tools, VLSI lithographic 
technology, and mathematical modeling of diffractive elements. 

Binary optics are constructed of micron-sized rectangular 
(high/low) relief patterns on dielectric or metallic surfaces. 
By controlling the depth, width, and period of the pattern, 
the amplitude and phase of an EM wave is transformed to 
produce the desired optical transfer function. The same tools 
used in VLSI electronic circuit fabrication are used to 
fabricate binary optics. 

The pattern formation is under computer control. Binary 
optics is named for the binary nature of making, etching, and 
coding technigues. Most of M.I.T.'s binary optics are made 
by contact printing the mask. The master lithographic mask 
is replicated onto a thin flexible glass mask, which is then 
vacuum contacted to the optical element so the flexible mask 
conforms to the surface, and then a contact print is made. 
However, for conventional optics aberration correction, the 
features are so large that one can simply projection print on 
curved or flat elements. Currently, M.I.T. Lincoln Labs is 
working on an extension of the technology which will allow the 
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diffractive patterns to be written directly onto the resist- 
coated optical elements with a laser beam, regardless of the 
curvature of the substrate. [Ref. 1] 

Binary optics are truly based on an existing 
infrastructure of fabrication technologies; it borrows from 
the VLSI electronics field. Additionally, binary optics 
nicely matches with the replication and embossing technology 
developed by the video and compact disc (CD) industry. New 
"plastics" and polymers are now transmissive in the IR over 
large bandwidths and could be used to mold whole systems, 
giving binary-optics technology low cost potential. [Ref. 1] 

A. HISTORY OF FABRICATION FOR BINARY OPTICS 

In the beginning (i.e., the 1960 's), a computer-generated 
hologram (CGH) , which is used as the mask for a binary optic, 
was plotted on a greatly enlarged scale and then photoreduced 
to obtain the working CGH. "The primary drawbacks to this 
method are plotter precision, resolution and speed, accuracy 
of the photoreduction steps, and the relatively long 
fabrication times due to the multistep process." [Ref. 5] 
Recently (i.e., the 1980's), electron beam (E-beam) 
lithography has been used to plot CGHs directly to scale, thus 
eliminating the photoreduction steps. The primary 
disadvantages of this method are: E-beam lithographic systems 
are very costly to buy and maintain, and the effective writing 
time of the CGH can be quite long. Since the first systems 
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developed were for printed circuit applications, they 
typically had a relatively large pixel size and writing area. 
Based on the success of the larger area writers, industry 
developed high resolution writers (i.e., laser writers to 
generate the mask) . [Ref. 5] 

Veldkamp has stressed that binary optics are not CGHs. 
Binary optics is a more comprehensive technology than CGHs. 
Computer-generated holography is concerned only with the mask 
design of optical elements but binary optics includes such 
design as well as highly efficient fabrication and coding 
technology. The production of highly efficient optical 
elements is what makes binary optics unique. In real systems 
applications, high efficiency elements are required. Veldkamp 
says "binary holograms" (i.e., a two-level phase hologram) 
are, at best, 41% efficient with each element causing a loss 
of 60% of the incident light, the net result being low signal 
levels at the detector due to poor efficiency. Consequently, 
binary phase holograms found few system's applications. In 
the Defense Advanced Research Projects Agency (DARPA) Binary 
Optics Program, Dr. Veldkamp and coworkers at M.I.T. Lincoln 
Labs learned how to make holographic relief elements highly 
efficient. [Ref. 1] 

High efficiency of binary optics is achieved by one of two 
methods. The first, relief structures can be constructed with 
a periodicity smaller than the wavelength of light but then 
the elements are very polarization and wavelength sensitive. 
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The preferred approach is to use structures larger than the 
wavelength of light, thus minimizing the polarization and 
wavelength effect. In the latter case, a staircase or multi- 
level pattern which is a repetitive mask application of the 
first method is used to achieve high diffraction efficiency. 
These staircase structures are used for broadband applications 
and M.I.T. Lincoln Labs have concentrated their research 
efforts in the last few years in this area. [Ref. 1] 

B. OVERVIEW OF FABRICATION PROCESS 

Figure 2.1 illustrates the overall design and fabrication 
process for making binary optics. 

First, the binary optic is designed using an optical ray- 
trace program; a very popular one used by industry is CODE V 
(Optical Research Associates) . CODE V generates a set of 
optimized phase coefficients defining the wave front phase map 
to be implemented, which is interpreted by in-house software 
for a particular pattern generator to define the mask for the 
lithography process. [Ref. 6] A primary constraint on the 
lens designer is that the minimum fringe spacing of the CGH 
does not exceed the resolution limit of lithographic mask 
generation hardware or alignment tolerances in the mask 
aligners. Additionally, the minimum period is equal to the 
minimum fringe spacing divided by the number of phase levels. 
Typically, E-beam machines are capable of eight micron fringe 
spacing. [Ref. 7] 

The making of a binary optic is similar to ordinary CGHs; 
binary optics use high-resolution lithography to encode the 



11 




Figure 2.1 Overall Design and Fabrication Process for 
Binary Optics [after Ref. 6] 
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fringe patterns defining the optical transformation performed 
by the element. The lithographic mask represents the CGH and 
typically yields diffraction efficiencies of around ten 
percent. However, binary optics use the mask only to 
delineate regions of the substrate that are to be etched after 
the element is properly processed with photoresist. A 
diffractive optic made from repeated mask applications (i.e., 
a staircase structure) has very high efficiency (>95%) . 

Three broad types of binary optics are illustrated Figure 

2 . 2 . 

Binary Elements 




1 - X h - X 



Fresnel Phase Lens 



\J 



T >> X 



X 



Mixed Optics 




Figure 2.2 Types of Binary Optics [from Ref. 6] 
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A binary optic may be a simple two-level structure, a 
multi-level structure, or a mixed optic (i.e., a binary optic 
with a conventional optic) . In general, binary elements can 
be made to operate for- minimum wavelengths in the IR, while 
Fresnel phase lenses (multi-level binary optics are made with 
repeated mask applications) can extend down to the visible; 
the limitations are from current lithographic machine 
resolution. In principle, the binary surface can perform 
arbitrarily large transformations on the incident wave front. 
However, resolution limits of current lithographic machines 
restrict the optical power (focusing ability) of the binary 
surface . 

1 . Mask Formation 

Various E-beam and optic lithographic mask-making 
machines are available with steps ranging from 0.1 to 1.1 
micron (jzm) ; a list is provided in Reference 6. 

a. E-beam 

Long exposures, in excess of 20 hours, can have 
unusually low yield factors. Perkin Elmer cites a yield of 
40% for a 25 hour exposure using their E-beam machine, 
Manufacturing Electron Beam Exposure System (MEBES-III) . A 
rough rule of thumb for the effect of pixel size on exposure 
time is that cutting the pixel size in half, may triple the 
exposure time. The exposure time lengthens since the number 
of stripes (lengthwise pass over the surface) increases by a 
factor of four and doubles the number of figures (shapes) 
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required to generate the pattern. The number of stripes 
covering the pattern to be exposed is directly related to the 
exposure time required, where stripe size is a function of 
pixel size. [Ref. 8] 

E-beam machines can use different geometrical 
figures in writing their patterns (i.e., parallelograms, 
rectangles, and trapezoids), and the choice of the shape 
influences the exposure time. Typically, rectangles are 
processed faster than parallelograms which are processed 
faster than trapezoids; the additional processing time is from 
the requirement to convert diagonal lines to a raster pattern 
while loading the stripe memory buffer. [Ref. 8] 

There are two major differences between the VLSI 
and CGH patterns: [Ref. 9] 

- The volume of data describing a CGH pattern is generally 
much greater than that of a VSLI pattern. 

- CGH patterns usually consist of a very small set of 
primitive shapes. 

Research labs and industry have developed their 
own in-house software to generate the mask from the optical 
ray-trace program's output to control the pattern writing 
machine. General rules to keep the size of the data to 
manageable proportions for the mask are: [Ref. 9] 

- The hologram pattern should generate the smallest number 
of primitive shapes (NPS) since the data is proportional 
to the NPS. Orientation of the patterns is also important 
(e.g., rectangular shapes oriented along the x-y axes 
generate less data than patterns consisting of curved 
lines in arbitrary directions) . 
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- The amount of graphical data is a function of the CAD 
application that generates the data. A tradeoff exists 
between processing speed and compact storage of data. 
Organization of small clusters of data into larger ones 
compacts the data, but lengthens the computer access time. 

- In VLSI applications (i.e., low data volumes), data 
fracturing programs are useful, but current CAD systems 
do not efficiently handle fracturing of CGH data. 
Therefore, special software is required for efficiency 
reasons . 

- Significant improvements in CGH design and fabrication are 
achieved if the CAD and E-beam machine are both available 
at the same location and connected together. CGH design 
and recording in parallel is thereby accomplished, thus 
reducing graphical data storage requirements and total 
fabrication times by significant amounts. Moreover, 
complex CGHs can be fabricated because the graphical data 
is generated in segments (i.e., as the E-beam machine is 
writing the previously generated segment, the CAD system 
can generate the next batch of data) . 

b. Laser writers 

Today, the majority of the patterns are generated 
using E-beam pattern generators. Many of the micro- 
applications require high resolution, requiring use of an E- 
beam generator. However, mask sets for 0.1 to 2.0 in binary 
diffractive optics can be rapidly generated on high resolution 
laser writers such as Texas Instrument's Laserad system. 
Currently, Texas Instruments has demonstrated the capability 
to use their laser writer to develop a simple Fresnel zone 
plate lens designed for use at 10 /xm wavelength. The 
advantage of using a laser writer compared to an E-beam 
pattern generator are: (Ref. 5] 

- The laser writer is less costly to operate and maintain. 

- The laser writer has a higher effective throughput. 
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Effective pixel sizes range from 10 unt to 0.5 fim 
for laser writers. Texas Instruments developed the Laserad 
system to plot photomasks for typical semiconductor 
applications. The Laserad I has position and overlay 
accuracies of better than 0.25 /zm, which is useful for 
generation of multi-level binary optics. 

2 . Reactive Ion Etching 

The fabrication of an N-level (multi-level) profile 
requires precision photolithography and processing. The masks 
for each relief level are used to contact print a similar 
relief profile in photoresist spun on the surface of the 
substrate. The processed photoresist consists of areas of 
bare surface or photoresist (unexposed areas) . After 
subsequent material removal by reactive ion etching, a profile 
of etched and unetched areas is produced. The photoresist 
acts as an etching barrier. Reactive ion etching allows a 
precise way of removing material anisotropically at the atomic 
level. In reactive ion etching, a plasma is created from the 
reactive atoms over the surface to be etched, thus enhancing 
the efficiency and precision of material removal. The process 
of producing photoresist patterns and then etching relief 
structures is repeated N times with N different masks, 
requiring complete registration, to produce an N-level relief 
(multi-level or staircase) pattern of different depths and 
heights. The key steps in fabricating a binary element are: 
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- Accurate registration of the N-level masks. 

- Precise control of depth and profile of material removed. 

Hence, the mask aligner and reactive ion etches are 
the critical pieces of equipment. [Ref. 7] 

3 . Embossing a Binary Optic in Plastic 

The embossing technique involves the following major 

steps: 

- The CGH is recorded in photoresist. 

- A master metal replica is formed by electroplating the 
photoresist pattern. 

- The metal master or a metal replica is used to repeatedly 
emboss its pattern into plastic. [Ref. 10] 

Polaroid fabricates the metal master from nickel. The 
developed photoresist plate is coated with a vacuum evaporated 
layer of silver, making it electrically conductive, and is 
immersed in an electroplating bath consisting of an 
electrolyte and a nickel anode. The resist plate acts as the 
cathode and a layer of nickel is plated onto the resist 
surface. Afterwards, the resist plate is stripped away, the 
nickel plate can be used for embossing. The nickel plate has 
a surface relief profile which is the mirror image of the 
original CGH pattern in the photoresist, making the pattern 
which is impressed into the plastic the same as the 
photoresist's. The plate can be mounted into a machine which 
allows for continuous embossing into plastic. The plastic 
must be softened by either heat, pressure, solvent or a 
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combination of these. Typical plastics used for embossing are 
vinyl, polycarbonate, mylar, or cellulose esters. [Ref. 10] 
Currently, M.I.T. Lincoln Labs, Polaroid, and 3M have 
demonstrated the ability to fabricate a binary optic in 
plastic. 3M Corporation mass produces commercially a binary 
optic used in a diffractive bifocal intraocular lens. [Ref. 
11] The replication (embossing) technology is the same used 
to manufacture CDs. Additionally new plastics transmissive 
in the IR are also available. [Ref. 1] If the lenses are 
produced in quantity, this technology is critical to 
lightweight and low-cost optics. 
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III. MILITARY AND SPACE APPLICATIONS 
FOR BINARY OPTICS 

The applications of binary optics to military and space 
systems is in its infancy; most current applications are still 
in the research and development stage. In the future, binary 
optics may be routinely used for; solar concentrators for 
solar panels, combining two-dimensional arrays of 
semiconductor lasers to produce a single powerful coherent 
laser beam, beam sampling, reshaping of laser beams for 
particular energy distribution, laser beam multiplexing for 
communication systems, beam steering, optical storage devices, 
head-up displays for pilots and astronauts, telescopes, and 
numerous other applications. 

Veldkamp divides binary optics into two classes of 
applications. Currently, the most attention is directed 
toward basic aberration correction with large features on flat 
and curved elements, and the direction is toward larger 
optics. The goal is to perform conventional aberration 
corrections: better, cheaper, and lighter. The second class 
of binary optics is comprised of hundreds of thousands of 
micro-optic elements of a few-inch-diameter wafer (e.g., 
microfine structures for delicate switching functions, 
multiplexers, filters, and amacronic sensors) . [Ref. 1] 
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This chapter summarizes current applications and 
speculates on the future for binary optics. Currently, both 
industry and research labs are trying to determine how to best 
exploit this new technology which has been called an "optical 
revolution," which will change the way future optical systems 
are designed. 

A. USING CONVENTIONAL OPTICS WITH BINARY OPTICS 

1. Equivalent Asphere 

A low-cost, conventional, spherical element can have 
a pattern etched onto the back surface and turn the lens into 
an equivalent asphere. Cost and producibility motivate this 
mix of conventional optics with binary optics. [Ref. 1] 

2 . Deliberate Dispersion 

Diffractive rings can be etched onto a surface in 
order to achieve a deliberate amount of dispersion into the 
element. Since this dispersion is opposite to most refractive 
dispersion, an achromatized element (i.e., all colors in a 
waveband focus at the same spot) can be made by balancing the 
mix of refraction and dispersion. [Ref. 1] 

Therefore, binary etch patterns allow control of both 
the surface profile of an element (asphericity) and the bulk 
properties of the material (dispersion) . 

3 . Correction of Spherical and Chromatic Aberrations 
When the staircase patterns are used to correct for 

spherical and chromatic aberrations in conventional elements, 
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the spherical curved surfaces of the conventional element 
performs the raw focusing power, and the features of the 
corrective diffraction patterns can be very large (e.g., 
typically hundreds of microns in width) . Since the features 
are larger than the wavelength of light, they become 
polarization-independent and color insensitive. Typical IR 
aberration corrections can be several hundred microns or more 
in width, and two or three microns in depth. However, 
depending on the application, features can also be submicron 
in scale. For example, fast planar optics (f/1) used as 
micro-lenses in front of detector arrays have the smallest 
lithographic feature as: 

_ 

S “ N" (3-D 



where 

s = size of lithographic feature, 

A = wavelength, 
f = f number of the lens, and 
N = number of phase levels. 

Therefore, features can be submicron in scale, depending on 
the application. [Ref. 1] 

4 . IR Applications 

Binary optics can improve the performance of IR 
systems when combined with refractive elements. Additionally, 
arbitrary phase profiles are possible since the elements are 
computer generated. Possible applications include forward 
looking infrared ( FLIR) systems, laser radar sensors, and 
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other laser applications (i.e., beam profile shaping and beam 
combining from diode arrays) . 

In general, optical systems are confined to spherical 
surfaces. Aspheric surfaces can be used to reduce the number 
of elements but are significantly more expensive to 
manufacture. The majority of IR lens elements are made from 
germanium, silicon, and zinc selenide. Unfortunately, these 
materials are inherently dispersive and therefore require 
optical systems with more elements to correct for the 
dispersion. However, by etching a diffractive surface onto 
a conventional refractive surface, the diffractive surface can 
correct for the dispersion and thereby reduce the number of 
required elements. The reasons for making a hybrid 
(combination diffractive and refractive) lens are: [Ref. 12] 

- The easily manufactured spherical surface can do the 
majority of the focusing. 

- The tolerances of making and aligning the diffractive 
masks are significantly reduced. 

- Spherical surfaces are inherently better in minimizing 
off-axis aberrations over a finite field of view. 

- The operative wavelength band of the system is greatly 
increased over that of a completely diffractive element. 

a. Broadband IR Imaging 

Information in this section is summarized from 
research by Fritz and Cox. [Ref. 13] 

Fritz and Cox investigated two optical designs for 
IR imagers, a fast staring system using an f/1 lens (field of 
6 degrees) and a slower scanning system using an f/2.4 lens 
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(field of 1.5 degrees). Their designs included hybrid 
refractive/diffractive elements as replacements for 
conventional spherical elements for use in the eight to twelve 
micron band. The diffractive element is used to apply 
aspheric or chromatic corrections. 

Mosaic arrays of uncooled detectors for low-cost 
staring systems require fast (<f/l) optics. Fritz and Cox 
found that the sensitive dependence of spherochromatic 
aberrations (introduced by the binary optic) on f number limit 
the performance improvements for a fast system (e.g., an 
uncooled staring system) . A f/1 triplet of conventional 
spherical elements cannot be replaced by a hybrid system of 
two spherical elements, each with an aspheric diffractive 
surface while maintaining resolution over the desired field 
of view and bandwidth. However, a 33% reduction in the number 
of elements could be achieved in the slower system using the 
hybrid elements with no loss in resolution, field of view, or 
bandwidth. 

B. IR FILTERS 

Dale Byrne has fabricated diffractive IR filters using 
microlithographic techniques; the minimum feature sizes less 
than 0.25 /xm for use as mesh filters with near and middle IR. 
These meshes may be used in applications like optically 
pumping far IR lasers which use a hybrid metal mesh dielectric 
mirror to have a high reflectivity at the pump's wavelength, 
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partial reflectivity and low absorptivity at the lasing 
wavelength, and couples a broad area of the laser mode. [Ref. 
14] Previously, three different mesh filters were used for 
ground based millimeter wave astronomical measurements and 
metal meshes were used as spectral filters in a far IR and 
submillimeter wavelength satellite radiometer. [Ref. 15] 

C. COHERENT LASER ARRAY BEAM ADDITION 

This section summarizes research by M.I.T. Lincoln Labs 
concerning coherent laser array beam addition. [Ref. 16] 

Coherent laser beam addition offers the potential to 
combine the power from many individual lasers into one 
extremely powerful coherent beam. The three techniques 
involved each use binary optics and can be employed to produce 
diffraction-limited laser sources from large two-dimensional 
arrays of laser diodes or by combining the output of several 
gas lasers. Specifically, the three techniques are: 
superposition using binary phase gratings, aperture filling 
by amplitude-to-phase conversion, and aperture filling by 
coupled microcavities. 

Currently, an individual diode laser is limited to 
approximately 25 milliwatts (mW) of power to prevent facet 
damage and device heating. Arrays of laser diodes overcome 
the heating problem by designing adequate space between lasers 
with the proper heat sinks; naturally, the power of each laser 
is limited in order to prevent facet damage. The collective 



25 



output of today's diode arrays is several hundred watts of 
quasi-continuous wave (CW) power with larger CW and pulsed 
arrays under laboratory development. Two-dimensional laser 
arrays have the potential of achieving power densities of over 
1 kilowatt per square centimeter. 

Power is extremely important if the laser is to be used 
as an illuminator, multimode fiber source, or laser pump. In 
some applications, power concentration and the far-field 
behavior is equally important; these applications include 
laser radar, optical communications, single-mode fiber 
sources, and optical printing. Arrays of mutually incoherent 
lasers have a low degree of spatial coherence across the 
array; therefore, if laser diodes are to be used for these 
applications, the outputs of the individual lasers must be 
combined coherently in order to achieve high luminance 
(luminance is a measure of power per unit area per unit solid 
angle and is sometimes called brightness) . To achieve high 
luminance and optimal far-field performance from laser diode 
arrays, three conditions must be met: 

- Mutual coherence must be established across the array 
insuring good spatial coherence for good collimation and 
focusing properties. 

- The phases of the lasers must be adjusted to provide the 
maximum power along the optical axis. Improper phasing 
results in beam broadening and reduced power along the 
optical axis. 

- The wave front emanating from the array is modified to 
produce a single lobed far-field pattern. The far-field 
pattern of an in-phase array consists of several narrow 
beams called beam combining. 
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It is of special interest to note that there are no 
conventional optic counterparts for the following methods of 
beam combining (i.e., binary optics are critical to the 
design) . 

1 . Superposition using Binary Phase Gratings 

The laser beams in the array may be combined via a 
binary phase grating to superimpose N lasing apertures, making 
the array effectively appear as a single emitter with N times 
the power density of single laser diode. The power per unit 
area increases by a factor of N while the divergence remains 
constant and the luminance increases by a factor of N, also. 
The laser beams are made to cross each other at a particular 
point in space, the resulting interference pattern is 
I E (x,y) | exp[ j (x,y) ) which is dependent on the beam angles and 
respective phases of the lasers. The binary phase grating is 
positioned in the interference plane to convert the field into 
a single beam. For perfect plane wave beam conversion, the 
transmittance of the grating, t(x,y), is the reciprocal of the 
interference, namely: 

4 x >y) = f E(x,y)| ex p{-M x ’y>} (3,2) 

The transmittance is less than or equal to one; any 
absorption will produce a loss in laser power, so the 
amplitude transmittance is chosen to be unity, and the 
amplitude component of the interference pattern is not 
corrected. However, the phase term is corrected by 
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constructing a grating to give the desired phase from Equation 
3.2. After the light passes through the conjugate phase 
grating, the light amplitude coupled into the zero order beam, 
A 0 , is given by the average of the resultant field: 



a ° ■ //: |E(x,y)|dxdy 



(3.3) 



The grating coupling efficiency is maximized by choosing 
the phases of the lasers that produce an interference 
amplitude |E(x,y) | which maximizes Equation 3.3. For example, 
six lasers are maximized using laser phases of n , n, 0, 0, n 
and tt , resulting in a coupling efficiency of 81% using 
Equation 3.3. 

The binary phase grating is placed inside a common cavity 
for the lasers (Figure 3.1) and combines the output into a 
single beam. Of particular interest is the fact that it can 




Figure 3.1 Superposition of Laser Diodes using a 
Binary Phase Grating [from Ref. 12] 
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be shown that the combining efficiency is identical to the 
efficiency of splitting a single beam into N separate beams 
(i.e., a good beam splitter is also a good combiner). 
Veldkamp used iterative techniques to optimize the desired 
diffraction pattern in a least squares sense for the 
combination of six laser beams; the six lasers were positioned 
along the +/- l, 2, and 3 diffraction orders of the grating. 

The angular plane wave spectrum of the grating output 
is shown in Figure 3.2 and the far-field output of the six 
lasers is compared with a single laser (scaled in intensity) 
in Figure 3.3. 

Although this technique has also been applied to HeNe 
lasers and C0 2 lasers, it is most appropriate for low 
divergence sources where additional collimation of the laser 
beams is unnecessary. High divergence sources (e.g., 
semiconductor lasers) require collimating optics; therefore, 
this method is limited to small arrays of semiconductor 
lasers . 

2 . Aperture Filling by Amplitude-to-Phase Conversion 

Aperture filling combines the laser beams by spreading 
the light uniformly across the entire laser array. The 
desired light distribution is uniform in amplitude and phase; 
the distribution is approximated by using binary elements to 
perform the required spatial filtering and phase correction. 
This method differs from superposition in that coherence is 
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Figure 3.2 Angular Plane Wave Spectra of Combined 
Laser Output [from Ref. 12] 
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Figure 3.3 Far-Field Output of Combined Beam Compared 
with Single Beam [from Ref. 12] 



30 



provided by the design of the array and not an external 
cavity. The array radiates in-phase so the phase component 
of the electric field is constant. However, the amplitude 
component of the fields are nonuniform due to the laser 
separation and their Gaussian mode structure; this amplitude 
variation produces unwanted grating lobes in the far-field. 
It can be shown that the percentage of energy in the central 
grating lobe equals the ratio of the area of the laser 
apertures to the area of the entire array; this ratio is 
defined as the fill factor. The remaining energy resides in 
higher order grating lobes. 

In this technique, two binary optical elements are 
required. The first element transforms the field by altering 
the phase to get the amplitude approximately constant and the 
second element, a phase-correcting plate, removes the phase 
variation introduced by the first. Both elements affect only 
the phase so no light is lost in the combining process; in 
theory, 100% of the light from an ideal array can be coupled 
into the central lobe. Computer simulations indicate coupling 
efficiencies of greater than 95% are possible with many real 
laser arrays. 

Figure 3.4 shows results from combining a ten-element 
AlGaAs array containing 2.4 micrometer laser sources spaced 
by six micrometers. Figure 3.5 shows aperture filling almost 
doubled the intensity in the main lobe (91% vice 51% of the 
total power produced by the array) . 
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Figure 3.4 Before Aperture Filling [from Ref. 12] 




Figure 3.5 After Aperture Filling [from Ref. 12] 
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The aperture filling technique yields a 
diffraction-limited far-field beam which is inversely 
proportional to the transverse dimension of the laser array 
vice the individual laser. Consequently, smaller diffraction 
spots and lower beam divergence is achievable without the use 
of beam-expanding optical elements; this feature allows 
simplification of some optical systems. This technique is 
limited to lasers with fill factors greater than 25%. 

Aperture filling may also be applied to beam combining 
of phased arrays of lasers for beam steering applications. 
Additionally, this method can be cascaded in order to transfer 
100% of the power into another system (e.g., use as a laser 
pump) . 

3 . Aperture Filling by Coupled Microcavities 

Diffractive microlens technology and Fresnel 
diffraction are used to create an array of small coupled laser 
cavities. The coupling locks the laser cavities together 
coherently and produces a narrow far-field pattern. The 
microlenses uniformly spread out the light via diffraction 
across the array, thereby producing the desired light into the 
central lobe of the far-field. The small apertures of 
semiconductor lasers have beams which expand into a cone of 
20 to 30 degrees and intersect with adjacent lasers a few 
hundred microns from the aperture. Collimating each laser 
beam at this point results in the desired uniform phase and 
quasi-uniform intensity of the wave front. 
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An array of diffractive optical elements is used which 
have the flexibility and precision required to adapt to 
virtually any size and shape. The lens requirements are quite 
strict but are achievable with binary optics, specifically: 

- The lens size must equal spacing of the lasers; typically 
ranging from ten to several hundred microns. 

- Each lens must be located as close to its neighbor as 
possible . 

- The shape of the lens aperture may be noncircular. 

- The lens should have low f number for efficient light 
collection . 

- Diffractive and absorptive losses must be low. 

- Each lens must produce an aberration-free collimated beam. 

- Each lens needs to have the ability to correct for 
aberrations in the original laser beam. 

a. Coupling Microcavities by Talbot Self Imaging 

Talbot discovered that illuminating a periodic 
object with coherent light produces an image at regular 
intervals due to free space diffraction. [Ref. 17] The 
diffraction pattern of the array has self imaging properties 
but each element also has diffracted light which continues to 
expand contributing to several elements in the first self 
image. This property is combined with Talbot imaging of the 
entire array to form a Talbot cavity. Light from a single 
laser is reflected from an external mirror which provides 
feedback to all the other lasers in order to ensure coherence. 
The output mirror is positioned so that one cavity trip 
corresponds to the Talbot distance, ensuring the free space 
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diffraction from the entire laser array produces a self image 
of the original laser apertures, allowing light to couple 
efficiently into the individual laser waveguides. If the 
lasers are not mutually coherent, the Talbot effect does not 
apply and no self image is formed. Consequently, much of the 
light falls between the laser apertures and is lost. Since 
the best efficiency only occurs for coherent operation, the 
system acts in a manner to sustain the coherent output. 

Lincoln Labs has used a linear array of AlGaAs 
diodes which used a Talbot cavity (Figure 3.6). The lasers 
are collimated using a microlens array and a Talbot cavity was 
formed by placing a planar mirror half the Talbot distance 
from the microlens array. 

The far-field pattern from six laser diodes is 
shown in Figure 3.7. The unlocked phase illustrates the 
output without the mirror positioned to satisfy the Talbot 
condition and the in-phase output illustrates the output for 
the mirror position satisfying the Talbot condition. The far- 
field pattern of Figure 3.7 was almost diffraction limited; 
85% of the energy was in the central lobe. Spectral analysis 
of the lasers showed lasing occurred for a common wavelength 
and all the lasers were locked in-phase. 

The beam shaping and aberration correction 
abilities of microlenses are particularly suited to combining 
large arrays of two dimensional diode lasers. Additionally, 
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Figure 3.6 Coupling Microcavity using Talbot 
Self Imaging [from Ref. 12] 
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Figure 3.7 Far-Field Pattern for Microcavity 
Coupling [from Ref. 12] 
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the entire system can be enclosed in a single substrate 
(Figure 3.8). The microlens array is etched into the front 
surface of the substrate, and a partially reflecting mirror 
is placed on the back surface. The thickness of the substrate 
satisfies the Talbot self-imaging condition. 



Figure 3.8 Microlens Array in Single Substrate 

[from Ref. 12] 

4 . Advantages of Laser Beam Addition 

Reliability is a key issue for any military or space 
system. Each of the methods of beam combination is immune to 
catastrophic failure since a single laser failure only 
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slightly diminishes the system's power output. The high 
efficiency and reliability of these applications to laser 
diode arrays makes them ideally suited to military and space 
applications . 

D. LASER BEAM PROFILE SHAPING 

This section summarizes research by M.I.T. Lincoln Labs 
concerning laser beam profile shaping. [Ref. 18] 

Veldkamp at M.I.T. Lincoln Labs has shaped a far-field 
laser beam profile from a Gaussian into a uniform intensity 
profile in either one or two dimensions. Uniform intensity 
profiles are needed in particular applications, such as laser 
radar systems utilizing heterodyning detector arrays, laser 
fusion, optical data processing, industrial laser heat 
treating of material surfaces, semiconductor annealing of 
material defects, vaporization of deposition materials, and 
improvements in energy extraction efficiency from high energy 
lasers. [Ref. 19] 

For power efficiency reasons, energy redistribution is 
more desirable than selective attenuation of the beam to 
achieve the desired profile. In a laser radar, the far-field 
intensity pattern must match the far-field projected image of 
the linear detector array, known as the radar footprint, with 
a uniform intensity. An N-fold beam compression was used to 
achieve the necessary far-field beam eccentricity for 
efficient coverage of the detector array's footprint. The 
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conversion of the Gaussian profile to a sine profile was 
achieved by spatial phase modulation of the beams using a 
binary diffraction grating. The grating produced constructive 
and destructive interference to help flatten the profile. 

The process is simply done in two steps. First, the laser 
beam intensity distribution is converted to a sine function. 
Next, the beam profile undergoes an N:1 anamorphic 
compression. Compression in the near-field results in 
expansion in the far-field in order to match the rectangular 
extent of the N-element array's footprint. As seen in Figure 
3.9, the shaper is composed of one phase grating and two one- 
dimensional anamorphic compressors. 

A far-field intensity profile for a C0 2 laser is shown in 
Figure 3.10. The dotted line shows the profile of the input 
beam, the output was scanned by a simple detector in steps of 
70 mils in the direction normal to the plane of Figure 3.10. 
Approximately 82% of the incident energy is contained within 
the first intensity crossings of the shaped intensity profile. 

E. LASER MULTIPLEXER 

Material in this section summarizes research by M.I.T. 
Lincoln Labs. [Ref. 20] 

Veldkamp has developed a binary diffraction grating which 
serves to multiplex local oscillator wave fronts for use in 
an IR laser radar system. The phase and amplitude matching 
of the signal to the local oscillator wave fronts maximizes 
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Figure 3.9 Diffractive Optics Performing Laser Beam 
Profile Shaping [from Ref. 13] 




Figure 3.10 Far-Field Intensity Profile of a TEMoo 
10.6 Micron C02 Laser with and without Beam Shaper 

[from Ref. 14] 
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the signal-to-noise ratio, provides the proper spatial 
resolution and minimizes the heat generated on the focal 
plane . 

The holographic surface relief grating takes a uniform 
plane wave from a single local oscillator (LO) laser and 
divides the beam into N plane waves focused on each of the N 
elements in phase with the incoming signal (Figure 3.11) . The 
virtual exit aperture calculated for the signal beam path is 
equal to the hologram size, resulting in a diffraction-limited 
LO amplitude distribution matched to the signal. As a result, 
the signal wave fronts and the LO wave fronts are spatial 
frequency bandwidth limited by the same low-pass filter, 
matching the signal precisely with the LO, resulting in 
diffraction-limited performance. Veldkamp believes it is 
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Figure 3.11 Holographic Generation of a Multibeam 
Local Oscillator [from Ref. 13] 
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feasible to illuminate a 100 element heterodyning detector 
array with LO beams from a single binary grating with 
reasonable power efficiency. 

This beam splitter performs a similar function as the one 
measured in Chapter IV, but the beam splitter in this section 
operates in the IR vice the visible waveband. 

F. SOLAR CELL CONCENTRATORS 

Currently, the Strategic Defense Initiative Organization 
( SDIO) is working on a survivable solar power array called 
Super and the Air Force is working on a Survivable 
Concentrator Photovoltaic Array (Scopa) ; both will eventually 
be used on military satellites. The Scopa designs include a 
solar concentrator; the purpose of the concentrator is to 
increase the power output of the cell by collecting more light 
and to increase the solar cell's survivability against enemy 
attack. Although binary optics are not currently considered 
for the concentrator's optics, there is no reason why they 
could not be used to perform the required function. 
Survivability against lasers is made possible since now the 
enemy must aim an antisatellite laser along the optical axis 
of the concentrator to get near the cell. [Ref. 21] 

The concentrating lenses focus sunlight at 25 to 100 times 
the normal intensity onto the cells. Concentrators are used 
because current and voltage in a cell increase linearly with 
the intensity of the sunlight that reaches them. 
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Concentrating the sunlight increases the power output of the 
cells and makes high efficiency possible. The use of 
concentrators results in higher temperatures in the solar 
cells, 75 to 100 C, which reduces efficiency. However, the 
gain in efficiency offsets the losses from increased 
temperature. Building cells to operate efficiently at those 
higher temperatures is still unrealized. Concentrators add 
weight, a critical factor in the design of space systems. The 
weight gain in the Scopa design is compensated by reducing the 
size of the solar cells; instead of using two by four 
centimeter cells, the new cells will be four to five 
millimeters in diameter. Additionally, fewer cells will be 
required. Another problem introduced by concentrators is the 
complication of designing the solar panels to track the Sun 
to within one degree to enable sunlight to enter along the 
optical axis of the concentrators; this problem remains to be 
explored. [Ref. 22] 

Polaroid has conducted research on a solar window or 
microscopic compound parabolic concentrator (CPC) . Polaroid's 
solar window is a holographic optical element (HOE) and not 
a binary optic; its fringe pattern is made by holographically 
interfering laser beams vice computer generating the pattern. 
It can be made transmissive selective with characteristics 
similar to a nonimaging CPC; it transmits the visible and near 
IR and reflects the longer wavelength IR. The CPC is a 
reflecting surface that is generated by rotating a parabola 
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around a central axis; it resembles a truncated paraboloid. 
Light enters the large end of the CPC over a fairly wide range 
of acceptance angles, is totally reflected within the CPC, and 
exits at the small end. The CPC can be made out of 
transparent dielectric material allowing total internal 
reflection at the walls in passing from one end to the other. 
The main advantage of the CPC is that it can track the sun 
over a large angle and still provide a reasonable amount of 
concentration; this may relax constraints on a solar panel 
control system for aiming the solar arrays. Polaroid has not 
considered this technology for space applications but there 
is no reason to prevent its application, but research for a 
suitable material is required. Polaroid has embossed 
microlenses of this shape into plastic from a master, 
demonstrating ease of manufacture and low cost. They 
constructed an array of lenslets with a diameter of 4 nm at 
the large end and tapering to 2 /urn at the small end (Figure 
3.12). If the space between the parabolas is either coated 
or filled with a highly reflective material, such as aluminum, 
then any long wavelength IR radiation that is incident on the 
embossed side will be reflected because the 2 /Ltm opening is 
too small for the radiation to pass. [Ref. 23] 

For space applications, further research is required to 
find an appropriate material suitable for use in the space 
environment. [Ref. 24] The space station will use deployable 
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Figure 3.12 Isometric View of Hexagonal Pattern 

[from Ref. 25] 

arrays of solar cells; in the future, a roll or collapsible 
panel of solar concentrators could be made and attached to the 
outside of such solar panels to increase the efficiency. 
Polaroid's ability to mass produce binary optics is an 
extremely important achievement because it will help 
demonstrate the low cost as well as the next generation of 
lightweight optics. 

G. CURRENT INDUSTRIAL APPLICATIONS 

Honeywell is currently researching the application of 
binary optics to ring laser gyro (RLG) readout optics, 
uncooled thermal imagers, helmet-mounted displays, 
antiref lective avionics displays, and optical storage disk 
readout. The RLG readout optics are for use in determining 
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the interference pattern from the missile's RLG, providing 25% 
cost reduction, a two fold improvement in signal to noise 
ratio, and easier assembly and alignment. The uncooled 
thermal imagers are using binary microlenses which provide a 
two-fold increase in sensitivity while reducing weight by 40%; 
they will be used on missile seekers and infantry night 
scopes. Helmet-mounted displays require lightweight optics? 
by using binary optics, a 25% reduction in weight, a 30% 
increase in the aperture, and reduced assembly and alignment 
costs were realized. [Ref. 26] The helmet-mounted displays 
currently designed are for helicopter pilots. The 
antiref lective avionics displays use a low cost plastic 
replication of a binary optic (estimated cost savings of a 
factor of ten) and provides reduced glare thereby enhancing 
safety. The optical storage disk readout may be used on 
aircraft as well as the space shuttle and space station. 
[Ref. 27] 

3M currently produces a diffractive bifocal intraocular 
lens which has been surgically implanted in patients in both 
the United States and Europe. The lens includes a blazed 
phase zone plate which directs most of the light to two 
diffractive orders. The combination of the diffractive and 
refractive surface enable the patient to have simultaneous 
bifocal vision, which is impossible with only conventional 
optics. The actual design is proprietary company information. 
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[Ref. 11] The diffractive lens is made using the replication 
procedure of embossing plastic from a master. [Ref. 28] 



H. FUTURE APPLICATIONS 

1. Near Term Future of Binary-Optics Technology 

Veldkamp speculated on the future of binary-optics 

technology in the next three to five years by saying: 

Two years ago, when we felt that the technology was 
developed far enough to be practical, we expected that it 
would be accepted by the optics industry with open 
arms — it could solve a specific but narrow set of 
problems. Quite the opposite happened. We went to 
missile seeker optics and other systems people. Basically 
they said: are you crazy? These are systems that were 

designed ten years ago, we're not going to change the 
design why should we? There's no money in it. So what 
is slowing down the acceptance of the technology is really 
the inertia in captive DoD (Department of Defense) 
industries. If we could put in the same technology 
transfer effort in the commercial optics market, 
acceptance would go ten times faster. There's a 
commercial market now for bifocal interocular lenses, for 
binary optics in fiber communication, for deep UV (ultra- 
violet) systems, for landscape lenses in compact folded 
cameras, for flat 3-D (three dimensional) displays, and 
for micro-optics in CDs and VCRs (video cassette 
recorders) . It is the commercial market that will have 
binary optics in place before it will appear in DoD 
systems on a large scale. [Ref. 1] 

2 . Amacronics 



The following is summarized from research by M.I.T. 
Lincoln Labs concerning amacronics. [Ref. 1] 

Amacronics is a layered structure of processing 
electronics, binary micro-optics, and detector arrays, with 
applications in imaging systems with focal plane processing. 

Today's sensors use a computer system that follows a 
conventional optics front end, the raster scanner system, and 
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the electronic detector array. A large computer is required 
to handle the signal processing for the array's output which 
is taken in a serial manner. Similar to computer systems, all 
electro-optic systems suffer from a "bottleneck"; computer 
systems use a single central processing unit which slows down 
the overall system. In electro-optical systems, all the 
optical information goes through a detector array at the focal 
plane which is the bottleneck. A massive amount of 
information comes in and is serially carried away from the 
focal plane which is in a raw form not easily handled by 
digital computers. 

Amacronic structures take advantage of parallel 
processing of the array output. Since processing occurs at 
the array, a smaller computer can be used to perform any 
remaining signal processing. Currently, M.I.T. Lincoln Labs 
is developing layered structures of optics and electronics in 
a parallel form (a processing unit per pixel) , similar to the 
human retina which uses similar amacrine clustered processing 
layers. "Amacrine" is from the Greek a macros meaning short 
range. Hence, the idea is to couple dynamically clusters of 
detector arrays. Similar to the vision of diurnal insects, 
amacrine structures may be able to dynamically allow tradeoffs 
between sensitivity and resolution. Carver Mead of Caltech 
has built networks tuned to detect motion with the electronics 
distributed in a uniform way. [Ref. 1] With binary optics, 
it may be possible to build systems with peripheral vision 
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much more motion sensitive than an on-axis fovial view, or 
systems tuned for edge detection or noise reduction. Dynamic 
control of detector arrays is the goal of amacronics. 

Military and space systems should see increased use 
of binary optics because they can save weight and improve 
performance of systems. In general, military applications 
will find increased use of binary optics because of their 
versatility and ability to accomplish what conventional lenses 
cannot do. 

I. ADVANTAGES AND COST 

1. Advantages of Binary Optics 

Binary optics can perform tasks that conventional 
optics cannot. One of the major advantages of binary optics 
is its adaptability to an optical system; the optic can be 
made in an arbitrary shape and be designed to yield the 
desired wave front while allowing very compact sensor designs. 
Additionally, using diffractive elements may reduce the number 
of required elements by 25 to 33 percent. [Ref. 29] 

2 . The Cost Issue 

A common claim is that binary optics are low cost. 
However, the cost issue is yet to be fully proved. [Ref. 27] 
3M and Polaroid have demonstrated the ability to emboss lenses 
in plastic from a master; if the lenses are produced in 
quantity, then the embossing technology is critical to low- 
cost optics as well as lightweight optics. If only a single 
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lens is to be produced, it is not clear if there is a cost 
advantage to using binary optics. However, if the system 
requires micro-optics (i.e., ten millimeters (mm) or smaller) 
in mass quantities or if quantities of at least 100 are 
planned for an equivalent system using binary optics, then 
binary optics should prove more economical and in some cases 
yield better performance. [Ref. 30] 
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IV. NEAR-ANGLE SCATTERING MEASUREMENTS 



Near-angle scattering is the light scattered less than one 
degree from the main beam. In optical systems, near-angle 
scatter degrades the image contrast, especially near bright 
sources or dark objects against bright backgrounds. [Ref. 31] 
As mentioned in Chapter I, the scattering measurements give 
a qualitative indication of fabrication errors and manu- 
facturing limitations. Knowing the amount of scatter from an 
optic indicates how well the optic will perform in the system. 
In the next chapter, the scattering measurements are compared 
to a scalar theory which is developed and modeled an a 
personal computer. 

The near-angle scattering from the specular transmitted 
beams of a binary optic beam splitter (transmissive optic with 
no focusing power) designed and manufactured by M.I.T. Lincoln 
Labs was measured. The optic used had known manufacturing 
defects (i.e., it was etched to an improper depth), and time 
constraints precluded obtaining a better specimen. The binary 
optic was designed to be a helium neon (HeNe) laser 
multiplexer, dividing a laser beam into seven separate equal 
intensity coherent beams (Figure 4.1). Because the binary 
optic was etched too deeply, the seven beams were not equal 
in intensity. The main lobe had slightly lower intensity than 
the first, second, and third lobes. Theoretical limitations 
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Third Lobe 
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First Lobe 

Main Lobe 
First Lobe 
Second Lobe 
Third Lobe 

Figure 4.1 Beam Splitter made by M.I.T. Lincoln Labs 

prevent the manufacture of an optic with only seven orders, 
hence unwanted orders will always be present. [Ref. 32] 

A. NEAR-ANGLE SCATTER INSTRUMENT 

The instrument used to measure the near-angle scatter was 
designed by Douglas W. Ricks at the Naval Weapons Center at 
China Lake, California [Ref. 4]. A top view of the instrument 
is shown in Figure 4.2. The instrument is capable of 
measuring scattered light from one degree to a few ten 
thousandths of a degree from the main transmitted beam. 
Scatter measurements at extremely small angles are possible 
because the detector is scanned through the specular beam 
which can then be subtracted from the scatter measurements 
(i.e., the background level may be subtracted from the 
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Figure 4.2 Near-Angle Scatter Instrument 
(Top View) [Ref. 31] 



measurement) . Small angular measurements are possible since 
the mirror has a long focal length and a very small horizontal 
scan distance (in the plane of Figure 4.2) is used for the 
detector. 




where 

/? = angle from the specular beam, 
d = horizontal scan distance, and 
f m = focal length of the mirror. 



(4.1) 
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By substituting d = 0.0002 in and f m = 71.6 in into 

Equation 4.1, the angle /3 representing the best angular 
resolution of the system is: 



0.0002 in 
71.6 in 



2.8 |iradians 



(4.2) 



Equation 4.2 assumes the sample has no focusing power. 

Additionally, the instrument is capable of measuring power 
down to eight or nine orders of magnitude from the peak of the 
transmitted specular beam. This is made possible by using a 
moderately high power laser, a sensitive detector, by 
eliminating stray light, by reducing the odd order aberrations 
through symmetry, and by carefully controlling the effects of 
diffraction . 

Lens samples are placed in one of two positions depending 
on whether the sample has a focal length. Figure 4.2 shows 
that a sample with no focal length is placed in the collimated 
beam between the mirrors. The binary optic beam splitter 
which was measured had no focal length, therefore the 
configuration of Figure 4.2 was used. If the sample is a lens 
with a focal length then it is placed after the mirror as 
shown in Figure 4.3 and the detector assembly is located at 
the combined effective focal length for the mirror and the 
sample. (The mirror shortens the effective focal length of 
the sample . ) 

A Spectra-Physics Stabilite model 124A HeNe laser, 
operating at 632.8 nanometer (nm) was the source. Following 
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Figure 4.3 Near-Angle Scatter Measurement for Sample 
with Focal Length (Top View) 

the laser is a plano-convex lens (convex side facing laser) 
with a 200 mm focal length. The ratio of the focal length of 
the mirror to the focal length of the lens determines the 
expansion of the laser beam. A large aperture of 7 mm on a 
15 by 15 centimeter (cm) plate was used following the 200 mm 
lens in order to eliminate stray light between the laser and 
the first mirror which is reflected by the optical bench and 
other surfaces. A small aperture, 500 /im, is placed at the 
common focus of the mirror and the plano-convex lens. The 
pinhole is approximately three or four times larger than the 
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focused beam at the 1/e points to avoid diffraction effects. 
(Each mirror and lens is at least four times the beam diameter 
in this apparatus.) The spherical mirrors are twins, each has 
a focal length of 71.6 in (1.82 m) and is approximately 12 in 
in diameter. 

The first mirror is placed to give the optimum collimation 
in the horizontal plane. The angle of incidence for the laser 
beam on each mirror must be the same in order to cancel the 
coma introduced by the spherical mirrors. Degradation occurs 
from spherical aberrations but is minimized by the long focal 
length of the plano-convex lens and mirrors. Since the 
spherical mirrors are used off-axis, astigmatism is introduced 
in the vertical plane. The astigmatism does not affect the 
system's horizontal resolution (i.e., the direction of scan 
is unaffected by astigmatism) . Although astigmatism could be 
completely eliminated by the use of off-axis paraboloid 
mirrors, they are costly and difficult to align. Normally, 
if the sample has no focal length it is positioned in the 
collimated beam between the two mirrors. A large aperture of 
7.5 cm on a 50 by 50 cm plate was used after the sample to 
block all stray light between the two mirrors. The binary 
optic holder actually blocks stray light as close as one half 
inch from the center of the binary optic. 

A pinhole is placed at the best focus in the horizontal 
plane of the second mirror. The silicon detector, made by 
Newport Corporation, is connected to their laser pico-watt 



56 



digital power meter, model 835. A 633 nm narrow-band filter 
is used in front of the detector; the filter has a 
transmittance of 57% at 632.8 nm and is used to block out 
ambient light. Translation stages with 0.0002 in resolution 
are used to position the pinhole and detector assembly at the 
focus. The pinhole and detector are simultaneously scanned 
in the horizontal plane to measure the scattered light from 
the sample. [Ref. 31] 

Laboratory air currents and temperature gradients can 
cause angular beam deflections comparable to the resolution 
of the instrument, so most of the instrument was enclosed in 
a Plexiglas enclosure. [Ref. 33] Two small windows were cut 
into the enclosure to permit the laser beam to enter unimpeded 
and to allow the experimenter to operate the micrometers on 
the translation stage. The laser was located outside the 
enclosure in order to isolate the heat generated from the 
laser with the instrument. 

The size of the pinhole in front of the detector involves 
a tradeoff between angular resolution and intensity available 
for the measurement. In other words, a small pinhole has 
excellent angular resolution but does not admit as much light 
so it has a lower signal-to-noise ratio than a larger pinhole. 
However, the small pinhole admits less stray scattered light 
and permits measurements at small angles. More than one size 
pinhole can be used to get a broader angular range of scatter 
information for a particular sample. 
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The instrument is optimized for an optic with no focusing 
power (i.e., the transmitted light is collimated). If the 
sample has a focal length, the configuration for the 
instrument in Figure 4.3 has less angular resolution. 
Additionally, the incoming beam to the sample is not 
collimated, and the focal plane becomes three dimensional 
which degrades results. A circular scan stage for the 
detector could compensate for this effect. 

B. PROCEDURE USED FOR MEASUREMENTS 

The actual beam splitter pattern is depicted in Figure 
4.4. From Figure 4.4, it is evident that the pattern is 
symmetrical; therefore, only four of the seven main orders 
were measured. 

The procedure to make measurements is relatively straight 
forward. First, a set of readings was taken without the 
sample in position using a particular size pinhole in front 
of the detector. These readings form the background readings 
which account for all the scattered light from the instrument 
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Figure 4.4 Actual Beam Splitter Pattern 
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and surroundings. Next, the sample was positioned and the 
measurements were repeated. The micrometer was moved the same 
distance relative to the respective peak; the angles are 
calculated by taking the arctan of the horizontal scan 
distance divided by the focal length. 

A lens with a focal length has less angular sensitivity 
since the effective focal length of the lens and mirror must 
now be used. The effective focal length is approximately 
equal to: 




where 

d = distance from mirror to lens, 
f 1 = focal length of the mirror, 
f 2 = focal length of lens sample, and 

f eff = effective focal length of mirror and lens combination. 
The effective focal length is always less than the 
focal length of the lens. 

C. SCATTERING MEASUREMENTS 

The data is presented in plots using a log-log scale 
because the scattered light varies over nine orders of 
magnitude. The lines on the graphs merely connect the data 
points . 

The abscissa has units of log sine angle of the scattered 
light in radians from the specular incident beam which is 
normal to the optic. Proceeding on the graphs from left to 
right, the angles become larger further away from the specular 
(nonscattered) beam. For very small angles, the sine of an 
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angle is the same as the tangent of that angle. Since near- 
angle scattering deals with angles less than one degree from 
the specular beam, the small angle approximation for the sine 
of the angle equalling the angle can be used. To illustrate 
the angular range of the instrument, consider the graph in 
Figure 4.5. The angular range is 10' 5 ' 6 < (3 < 10’ 2 ‘° radians 
which corresponds to 0.00014 < (3 < 0.573 degrees. Each beam 
starts with the specular beam as zero reference for the 
angular measurements. 



Background 




Figure 4.5 Background Measurements for 25 and 500 

/nm Pinholes 

The term Bidirectional Transmission Distribution Function 
(BTDF) was invented by Edwards and is universally used by 
professionals [Ref. 34]. The bidirectionality arises from the 
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incidence and exitance directions [Ref. 35]. Exitance is the 
radiant flux density of a wave field leaving a surface at a 
normal direction. The ordinate is in terms of the log of the 
BTDF, which has units of 1/steradian. The BTDF is: 

BTDF = ~ ^ al ^ (4.4) 
Pin 



where 

P scat (/3) = power scattered as a function of (3 per solid angle 
in steradians (sr) , 

/? = scattering angle from the specular beam, and 

P jn = power incident on the optic, which is constant. 

The solid angle n subtended by the detector as seen by the 

scattering test optic is given by: 



Q = 




(4.5) 



where 

A p = area of detector pinhole, and 

f m = focal length of mirror (for test optic with no 
focusing power) . 

If P det (/?) is the power scattered through the detector at 
angle then 






pjp) 

n 



Substituting into Equation 4.4, the BTDF is: 



(4.6) 



BTDF = 




(4.7) 



To illustrate the calculation of the BTDF, values from the 
experiment are substituted into Equation 4.7: 
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BTDF = 



(1.82 m) 2 
1 25 \im\ 



/0.26mW\ 
\ 6.3 mW J 



(4.8) 



BTDF = 2.8 * 10 8 (4.9) 

Log (BTDF) = 8.4 (4.10) 



The 25 /im and 500 /im pinholes were measured and verified 
to be the proper size using a differential interference 
contrast (Nomarski) microscope. 

To adjust data for displaying two different size pinholes 
on the same graph, the respective BTDFs are set equal to each 
other : 
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Isolating P 25 in Equation 4.11 yields: 



P25 = 



A25 

A 50 O- 



P500 



(4.12) 



where 

P 500 = P owe r readings for 500 ism pinhole, 

P 25 = power readings for 25 /im pinhole, 

Ajqo = area of 500 /im pinhole, and 
A 25 = area of 25 ism pinhole. 

The results show that P 500 is multiplied by A 25 /A 500 

1 . Background Measurements 

The background measurements without the sample present 

are shown in Figures 4.5 and 4.6. An expanded view of 10’ 3 ' 5 

to 10’ 2 radians is shown in Figure 4.6 where it is evident that 
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Figure 4.6 Expanded View of Background Measurements 

the larger pinhole tends to average the readings. The range 
of the instrument from Figure 4.5 is slightly more than nine 
orders of magnitude for the 25 /im pinhole and six orders of 
magnitude for the 500 /xm pinhole. The data for the larger 
pinhole starts at 10’ 3,6 radians because a larger pinhole lets 
in so much light that the detector saturated for an unimpeded 
specular beam (i.e., background measurements). The larger 
pinhole lacks angular resolution closer to the specular beam 
since the detector views a larger area and thereby tends to 
average light over a larger area. Two sizes of pinholes 
provide better information. The larger pinhole gives averaged 
reading further from the specular beam, and should have a 
larger signal-to-noise ratio for low scattering levels. 
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Background measurements are presented to indicate the 
angular range in which the background may be considered 
insignificant. If insignificant is defined as three orders 
of magnitude from the specular peak, then readings for angles 
greater than or equal to 10' 4 radians can be considered to have 
little influence by the background, whereas readings in the 
angular range from 10’ 5 to 10’ 4 radians would be influenced by 
the background. One may choose a less or more stringent 
definition for insignificant levels of background. Background 
measurements should be compared to measurements of the 
scattering test optic. If the relative background scatter is 
nearly equal to or greater than the relative scattering from 
the optic being tested for scatter, then the background 
scatter is significant. If the relative scattering at some 
angle from the optic under test is much greater than the 
relative background scattering at that angle then the 
background scattering is insignificant. There is no standard. 

2 . Main, First, Second, and Third Lobes 

The experimentally measured main, first, second, and 
third lobes are depicted in Figures 4.7 through 4.14. Figures 
4.8, 4.10, 4.12, and 4.14 are expanded views of the respective 
lobes at larger angles. Readings for each lobe start to 
increase at the larger angles since light is received from the 
adjacent beam (see Figure 4.4). Again, the 500 jum pinhole 
averages the readings and the same conditions noted with the 
background readings occur. The background and beams were 
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Figure 4.7 Main Lobe Measurements for 25 and 

500 /im Pinholes 
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Figure 4.8 Expanded View of Main Lobe 
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Figure 4.9 First Lobe Measurements for 25 and 

500 /im Pinholes 




Figure 4.10 Expanded View of First Lobe 
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Figure 4.11 Second Lobe Measurements for 25 and 

500 /xm Pinholes 
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Figure 4.12 Expanded View of Second Lobe 
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Figure 4.13 Third Lobe Measurements for 25 and 

500 Pinholes 
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Figure 4.14 Expanded View of Third Lobe 
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measured by moving the micrometer 0.0002 in for 100 data 
points, 0.001 in for 100 data points, and then 0.01 in for the 
remaining data points (approximately 25 to 100 additional 
points). The detector 'stage was moved in 0.005 in increments 
for the 500 /im pinhole instead of the smaller increments used 
previously (i.e., 0.0002, 0.001, and 0.01 in); approximately 
75 data points were obtained. Fewer data points are required 
for the larger pinhole since the readings do not change 
significantly with small horizontal detector movement, 
a. Main Lobe 

The measurements for the main beam in Figure 4.7 
show the range of the measurement as approximately six orders 
of magnitude; the range is smaller than the background due to 
the scattering which dominates at larger angles. The 
background shows a range of about eight or nine orders of 
magnitude, indicating approximately two orders of pure scatter 
in the main lobe. The power in the beam measured for the 
background is now divided into the separate lobes (see Figure 
4.4) and the specular power level is reduced. A small 
increase in the 25 /zm readings for angles greater than or 
equal to 10' 3 ‘ 75 radians is evident in Figure 4.7 (i.e., the 

slope does not have a gradual decrease as does the background 
in Figure 4.5); the irregularity in the slope is probably 
associated with a manufacturing defect. As mentioned earlier 
the main beam is of slightly lower intensity than the first, 
second, and third lobes which is attributed to the improper 
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etch depth. The upward trend at 10' 2 ' 5 radians is due to the 
received power from the adjacent first lobe (see Figure 4.4). 

b. First Lobe 

The measurements for the first lobe in Figure 4.9 
show the range of the measurement as approximately seven 
orders of magnitude. As before, greater resolution is 
obtained with the smaller pinhole and the larger pinhole 
averages the readings. The data closely follows the shape of 
the background reading indicating good performance. The 
background shows a range of about eight or nine orders of 
magnitude, indicating one or two orders of pure scatter in the 
first lobe. Power from the adjacent second lobe is detected 
at 10' 2,5 radians producing an upward trend. 

c. Second Lobe 

In Figure 4.11, the second lobe displays 
characteristics similar to the main lobe. Slightly less range 
is seen when compared to the first lobe, approximately 6.5 
orders of magnitude vice 7. The background shows a range of 
about eight or nine orders of magnitude, indicating 
approximately two orders of pure scatter in the second lobe. 
The intensity does not fall off as smoothly as the background 
(see Figure 4.5). Manufacturing defects are probably 
responsible for the increase in intensity after 10" 3 ‘ 5 radians. 
The upward trend at 10' 2 ’ 5 radians is from received power in 
the third lobe. 
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d. Third Lobe 



The third lobe measurement, Figure 4.13, displays 
similar behavior to the main and second lobes. The range is 
6.5 orders of magnitude for the 25 /nm pinhole. The background 
shows a range of about eight or nine orders of magnitude, 
indicating approximately two and a half orders of pure scatter 
in the third lobe. Again, an increase in intensity is 
observed for angles greater than 10’ 3 " 6 radians probably due to 
manufacturing defects. The intensity does not rise as much 
at larger angles since the next adjacent lobe is an unwanted 
lobe (see Figure 4.4). 
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V 



SCATTERING THEORY FOR BINARY OPTICS 



A scalar scattering theory for binary optics was 
originated by J. Merle Elson at the Naval Weapons Center, 
China Lake, California [Ref. 36]. The author collaborated 
with Elson in developing the theory? the author modeled the 
theory using a personal computer and compared the theory with 
the measurements from Chapter IV. This relatively simple 
theory was modeled using Microsoft FORTRAN and run on an IBM 
compatible personal computer. The scalar theory for binary 
optics yields approximate results; exact results require a 
full EM theory which considers boundary conditions and the 
polarization of the light. More accurate results are computer 
time-intensive, costing thousands of dollars to run a 
simulation. Additionally, the complexity of programming or 
running an existing program has increased by several orders 
of magnitude. For these reasons, a scalar scattering theory 
is a more attractive alternative. 

The scalar scattering theory developed considers 
scattering from two types of surface topography. Scattering 
is produced from surface roughness (micro-irregularities) and 
large scale surface features (i.e., binary step pattern). The 
micro-roughness is assumed to be generated by a random process 
and the binary surface profile is assumed to be deterministic. 
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The theoretical model developed is compared with 
experimental near-angle scattering measurements from Chapter 
IV. 



A. SURFACE ROUGHNESS MEASUREMENTS 

Surface roughness may be measured by using an instrument 
that gives a surface profile by either one of two methods. 
One method involves contact with the sample and the other does 
not. The contact method uses a probe which contacts the 
surface and follows height variations as either it or the 
surface is translated. The noncontact method uses optical 
interferometry techniques to measure surface height 
variations. Advantages and disadvantages exist for both 
methods and are thoroughly discussed in Reference 37. 

A contact method using a Talystep instrument was used to 
measure the roughness of the binary optic. The Talystep uses 
a diamond styles probe that touches the surface. Talystep 
measurements were performed by Jean Bennett at the Naval 
Weapons Center, China Lake, California. Figures 5.1 and 5.2 
show the Talystep profile of the binary optic measured. The 
average height of the optic is 6554 angstrom (A) and the 
period is 200 ^ m . A small defect is observable in one of the 
trailing edges of the binary optic and is more evident in 
Figure 5.2. The Talystep can measure roughness less than 1 
A root mean square (rms) on polished surfaces. [Ref. 37] The 
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horizontal resolution of the Talystep is determined by the 
diameter of the stylus (0.6 /xm for the measurements in this 
chapter) . 

Figures 5.3 and 5. -4 show the roughness of the top and 
bottom of a groove measured using the Talystep. In Figures 
5.3 and 5.4, 3.5 A roughness was measured on the top and 4 A 
roughness was measured on the bottom. The binary optic was 
made using reactive ion etching which produced an extremely 
smooth surface. 

The instrument noise for the Talystep was 0.5 A rms . 
Figure 5.1 used 13,098 points on a 500 profile with a 
0.0382 /xm sampling interval. The profile length for the 
groove top was 7 0 /xm and the bottom was 3 0 /xm. 

The autocovariance function is a measure of the correla- 
tion properties of the random surface roughness. It is the 
product of two "copies" of the same surface profile where one 
copy is shifted relative to the other. The amount of lateral 
shift between the two copies is called the lag length. A 
large (small) value of the autocovariance function indicates 
that the random surface roughness is statistically correlated 
over larger (smaller) lag lengths. In other words, a large 
value of the autocovariance function indicates that one "copy" 
is a good approximation for the other "copy" for that 
particular lag length. The value for zero lag length is 
important because it equals the square of the rms roughness 
of the profile. The correlation length may be defined as the 
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Figure 5.3 Groove Top Surface Roughness 




Figure 5.4 
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value at which the autocovariance function drops to 
1/e (0.368) of its value at zero lag. Alternate definitions 
use 1/10 of the value at zero lag, or where it first becomes 
zero. The autocorrelation function is sometimes used instead 
of the autocovariance function. The autocorrelation function 
is normalized to have a value of unity for zero lag length and 
is dimensionless. [Ref. 37] 

From the Talystep data, the autocovariance function was 
computed to verify the form of function assumed for the 
autocorrelation function. Since the autocovariance function 
is normalized to one, it is the same as the autocorrelation 
function. Figures 5.5 and 5.6 show the average autocovariance 
for the top and bottom of the grooves. 

An exponential autocorrelation function is a traditional 
assumption that appears to give good results in many cases. 
Unfortunately, the autocorrelation length from the Talystep 
can be misleading due to the short length of the scan run. 
Near-angle scatter is primarily concerned with long 
autocorrelation lengths on the surface; therefore, a long scan 
run for the Talystep is desirable. Previous measurements made 
by Wyco Corporation on this binary optic using a Wyco 
profilometer which uses optical interferometry techniques, 
yielded autocorrelation lengths ranging from 1 iim to 1,000 jum* 
[Ref. 38] 
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B. SCALAR SCATTERING THEORY 



The scalar theory developed assumes a Gaussian wave with 
constant phase incident on the binary optic. The phase 
changes after the aperture or source plane (see Figure 5.7) 
are due to the step height variations (i.e., the binary 
surface) and the surface roughness. It should be emphasized 
that this is a scalar solution and vector notation is used as 
a convenient shorthand in the following notation (i.e., the 
phase and polarization of the light are not considered in the 
theory presented) . 



A u(p.'t') 



Incident 

Wave Source Plane 

(Aperture) 




Source Plere 

■< 2 



U(p\ 4/) 



Observation 

Plane 




> 



Figure 5.7 Diagram for Development of Theory 
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In Figure 5.7, z is the distance between the source and 
observation plane, and U is the intensity which has both 
amplitude and phase. The field after the aperture assumes no 
absorption within the optic or back scatter from light 
reflecting off the incident face of the binary optic. The on- 
axis incident amplitude is assumed to be unity with a Gaussian 
amplitude distribution. 

For Fraunhofer diffraction the diffracted wave is: [Ref. 



39] 

u(x\y') 



I ik L ,2, ..,2j \ 



exp(ikz}exp^[x' +y' JJ 
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-f CO 

U(x,y) exp 

CO 



f-J2n 

\ Xz 



[xx'+yy'] 



dxdy 



(5.1) 



where 

U(x,y) = amplitude in the source plane, 

U(x',y') = amplitude in the observation plane, and 
z = distance from the exit pupil to the observation 

plane . 

In the experiment, the second spherical mirror focuses the 
collimated beam at the location of the detector which 
eliminates the quadratic phase term. 

To account for phase changes due to the binary step 
pattern, a term accounting for the phase change is added. 
Assuming additional aberrations are introduced by random 
surface roughness, £ (x,y) , an additional phase term is added 
to model this roughness. A Gaussian amplitude term is used 
for U(x,y) since a TEM 00 laser was used as the source. These 
modifications are incorporated into Equation 5.2. 
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Fourier series expansion of phase changes due 
to the binary height variations of surface, 

phase change due to the random roughness 
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= random roughness on surface, 

= index of refraction (constant) , 

= period of binary optic (constant) , 

= Fourier coefficients (calculated for binary 
surface 
profile) , 

= wavelength (constant) , 

= radius of the beam width (constant) , and 



= Gaussian laser beam with unit amplitude. 



The integral in Equation 5.2 is evaluated in Appendix A, 
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(5.6) 
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C. DETERMINATION OF FOURIER COEFFICIENTS 

To evaluate the Fourier coefficients, a p , the surface 
profile, h(x) , must be known. To determine h(x) for the 
binary optic, measurements were made and then averaged from 
a photograph of the optic and a ruler at the same 
magnification using a Nomarski microscope (Figures 5.8 and 
5.9). Each increment on the ruler in Figure 5.9 corresponds 
to 1 jLxm. The surface profile shown in Figure 5.10. 

The height, H, was previously determined using the 
Talystep as 6554 A. From the photographs, X, = 42.92 iim and 
X 2 = 8 6. 5 urn. 

To analytically find the Fourier coefficients, we use: 



h(x) = surface profile as a function of x, and 
a p = Fourier coefficients. 

The integral in Equation 5.7 is evaluated in Appendix B, 
resulting in: 

If p = 0, then 




(5.7) 



where 



n 



index of refraction 




(5.8) 
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Figure 5.8 Photograph of Binary Optic Using a 
Nomarski Microscope 




Figure 5.9 Photograph of Ruler Using a 
Nomarski Microscope 
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Figure 5.10 Surface Profile of Binary Optic 



If p f 0 f then 




Substituting values from the experiment into 5.8 and 5.9 
will give the Fourier coefficients. The period, L, is 200 /im; 



the height, H, is 6554 A; X, is 42.92 /im; X 2 is 


86.5 


/im ; and 


the index of refraction, n, for quartz is 


approximately 


1.45702. The square of the magnitude of 

coefficients is: 


the 


Fourier 


| a 0 | 2 = 0. 0219 




(5.10) 


| a 1 1 2 = 0.128 




(5.11) 


| a 2 | 2 = 0.141 




(5.12) 


| a 3 | 2 = 0.136 




(5.13) 
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More coefficients may be computed but they will have 
negligible effect on the calculation since these Fourier 
coefficients correspond to the measured beams of the binary 
optic and the magnitudes are later squared which makes their 
effect even smaller. 

D. COMPUTER MODELING 

The results from 5.3, 5.8, and 5.9 were programmed using 
Microsoft FORTRAN for an IBM personal computer. (The program 
requires a math co-processor.) A copy of the program is 
included in Appendix C. The data files generated by the 
program were plotted along with experimental values for the 
25 and 500 fj,m pinholes using Hewlett Packard's Charting 
Gallery program on an IBM personal computer. 

The following assumptions were used in modeling the 
theory. The form of the autocorrelation function is 
exponential (traditional form for function) , the average 
surface height is 6554 A, the rms surface roughness is 4 A, 
and the period is 200 fxra (where the average surface height and 
rms surface roughness are from the Talystep measurements, and 
the surface period was measured from photographs) . 

E. COMPARISON OF EXPERIMENT WITH THEORY 

The theoretical values generated by the computer program 
in Appendix A are compared with experimental data. Two sets 
of graphs were produced. First, a set of graphs was produced 
to yield a best fit of the theory for each lobe to the 
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experiment by adjusting values for autocorrelation length and 
roughness (Figures 5.11 through 5.14) . A second set of graphs 
was produced for each lobe using the same value of roughness 
and autocorrelation length (Figures 5.15 through 5.18). The 
lines on the graphs merely connect the data points with the 
exception of the solid line labeled "theory" which is obtained 
from theoretical equations. 

The data is presented in a log-log format because the data 
varies over nine orders of magnitude. The abscissa has units 
of log sine angle of the scattered light in radians from the 
specular incident beam which is normal to the optic. 
Proceeding on the graphs from left to right, the angles become 
larger further away from the specular beam. Since near-angle 
scattering deals with angles less than one degree from the 
specular beam, the small angle approximation for the sine of 
the angle equalling the angle may be used. To illustrate the 
angular range of the instrument, consider Figure 5.11. The 
angular range is 10" 3 ' 85 to 10' 2 ' 55 radians which corresponds to 
0.00809 to 0.16148 degrees. The ordinate is expressed in 
terms of log relative intensity; the data is normalized to the 
power in the specular peak for each beam. 

1 . Best Fit Graphs (Figures 5.11, 5.12. 5.13, and 5.14) 
a. Main Lobe 

In Figure 5.11, a best fit for the main lobe was 
obtained using a correlation length of 0.75 mm (750 jum) and 
roughness of 250 A. The correlation length from Figure 5.5 
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Figure 5,11 Comparison of Theory for Main Lobe 
(Autocorrelation Length 750 /zm, Roughness 250 A, Best Fit) 



25 500 Theory 
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Figure 5.12 Comparison of Theory for First Lobe 
(Autocorrelation Length 100 /zm, Roughness 50 A, Best Fit) 
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Figure 5.13 Comparison of Theory for Second Lobe 
(Autocorrelation Length 200 Jim, Roughness 7 5 A, Best Fit) 
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Figure 5.14 Comparison of Theory for Third Lobe 
(Autocorrelation Length 200 ^m, Roughness 125 A, Best Fit) 
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Figure 5.15 Comparison of Theory for Main Lobe 
(Autocorrelation Length 175 /im, Roughness 80 A, Average Fit) 
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Figure 5.16 Comparison of Theory for First Lobe 
(Autocorrelation Length 175 jxm. Roughness 80 A, Average Fit) 
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Figure 5.17 Comparison of Theory for Second Lobe 
(Autocorrelation Length 175 /xm, Roughness 80 A, Average Fit) 
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Figure 5.18 Comparison of Theory for Third Lobe 
(Autocorrelation Length 175 ;xm, Roughness 80 A, Average Fit) 
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is 1 /m at the 1/e point but previously measured correlation 
lengths from the Wyco profilometer ranged from 1 n m to 
1,000 Am so 750 /m is reasonable. [Ref. 37] The roughness 
was much higher than the Talystep indicated (4 A from Figure 
5.4) . 

_ V or 

The scatter in the angular range of 10 ' to 
10-3.68 S p ecu i ar scatter from the shape of the binary 
optic's surface and has a constant slope. It is evident that 
the theory predicts less scatter in this region. The scatter 

.7 /O , 

at angles larger than 10 is due largely to the micro- 
roughness on the surface; this area on the graph departs from 
the constant slope indicating the scatter due to roughness is 
overpowering the scatter due to the binary optic's shape. The 
theory agrees will with the experiment in the larger angular 
range. The upward trend of the data at 10’ 2 ' 6 radians is from 
the detector sensing energy from the next lobe. 

b. First Lobe 

In Figure 5.12, more reasonable values were used 
for the autocorrelation length and roughness, 100 nm and 50 
A respectively. The theory predicted less scatter in the 
angular region of 10' 3 ' 85 to 10’ 2 ' 5 radians. Overall, the theory 
is in good agreement with the experiment. 

c. Second Lobe 

A best fit for the data in Figure 5.13 used 200 
Aim correlation length and 75 A roughness. The theory fit the 
data well for the angular range of 10’ 3 ' 4 to 10~ 2 ' 55 radians. 
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d. Third Lobe 



A best fit for the data in Figure 5.14 used 200 
/xm for the correlation length and 125 A roughness. Similar 
fit to the second lobe was obtained. 

2 . Average Fit Graphs (Figures 5,1s, 5.16, 5.17, and 

5.18) 

Based on the best fit values for roughness and 
correlation length used to generate Figures 5.11 through 5.14, 
average values were obtained. An autocorrelation length of 
175 /urn and roughness of 80 A was used for each lobe in Figures 
5.15 through 5.18. 

a. Main Lobe 

In Figure 5.15, moderate discrepancies exist 
(i.e., an order of magnitude difference in the angular range 
Iq- 3.34 10 ' 3 ‘° radians) where the theory previously fit the 
data well. The improper etch depth in the binary optic caused 
the main lobe to be lower in intensity than the other lobes 
and may account for the variation not present in the other 
lobes . 

b. First Lobe 

Overall, a good fit of the data in Figure 5.16 is 
obtained, the theory closely follows the experiment from 
IQ' 3 ’ 45 to 10' 2 - 5 radians. 
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c. Second Lobe 

In Figure 5.17, a good fit of the data is 
obtained. Theoretical results coincide well with the 
experiment for 10" 3 ’ 4 to 10' 2 ' 55 radians. 

d. Third Lobe 

In Figure 5.18, results are similar to the main 
lobe, moderate discrepancies exist but the shape of the theory 
curve matches the shape of the experimental curves. 

F. DISCUSSION OF RESULTS 

The graphs in Figures 5.11 through 5.14 are best fits of 
the theory to the experiment by adjusting the autocorrelation 
length and surface roughness in the computer model . The 
average fit in Figures 5.15 through 5.18, which use the same 
value for the autocorrelation length and surface roughness, 
yielded poorer results, especially in the main lobe. Clearly, 
the values of autocorrelation length and roughness obtained 
from the Talystep measurements were either unrepresentative 
of the binary optic or the scalar theory requires further 
revision. The discrepancies in using larger values for the 
autocorrelation length and roughness may be attributable to 
the short scan length, 500 /xm, used by the Talystep, which is 
a physical limitation of the instrument. 

The short scan length used by the Talystep leads to two 
problems. The first problem concerns the instrument's 
inability to measure the overall roughness across the entire 
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optic. Therefore the Talystep may not be able to determine 
long autocorrelation lengths which primarily affect near-angle 
scattering. The Talystep is merely giving us a sample of the 
binary optic's surface, although runs are made in different 
areas of the optic. The question remains exactly how to 
determine the roughness and autocorrelation length. The 
variety in the Wyco's autocorrelation measurements are also 
frustrating. Perhaps the roughness and autocorrelation length 
vary with both the scan direction and the length of the run 
on the surface of the optic. The Wyco profilometer looks at 
different surface areas depending on the magnification of the 
lens used and will consequently generate different values. 
Additionally, the Wyco profilometer is highly dependent on the 
experience and skill of the operator, especially when using 
low magnifications which view a larger surface area of the 
optic . 

In general, the graphs indicate that the theory conforms 
well to the experiment. The scalar theory follows the shape 
of the actual measurements, especially in the central angular 
region. The greatest disagreement occurs at the smallest 
angles, where the measurements exceeded the theory by several 
orders of magnitude. 

The average fit for the graphs in Figures 5.15 through 
5.18 indicates some disagreement between the theory and 
experiment for the main and third lobes. The worst fit occurs 
for the main lobe which may be attributable to the improper 
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etch depth of the optic which caused the main beam to be lower 
in intensity. The longer autocorrelation lengths may be 
justified by using the Wyco profilometer instead of the 
Talystep. The Talystep was not desirable due to the short 
profile lengths. 

The results are very promising for using a scalar theory 
to predict the scattering from binary optics. The theory is 
relatively straight-forward and precludes time-intensive 
computer simulation using the full EM theory. The scalar 
theory is a good approximation but it did not account for: 

- Back scatter off the front (incident) surface of the 
binary optic (i.e. ; all incident light is transmitted). 

- Any inhomogeneities or defects present within the material 
the binary optic is made from. 

- Polarization of the light. 

Future work should determine how to best obtain 
measurements for the autocorrelation length and roughness of 
a binary optic. It may also be helpful to revise the 
theoretical model to include a random offset in the period of 
the binary step pattern. An attempt to account for the 
backscatter or polarization would complicate the model 
excessively. 
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VI 



SUMMARY 



A. ADVANTAGES AND APPLICATIONS FOR BINARY OPTICS IN MILITARY 
AND SPACE ELECTRO-OPTICAL SYSTEMS 

In summary, binary optics will never totally replace 
conventional optics. However, by using binary optics the 
following may be realized: 

- Correction of spherical and chromatic aberrations in 
conventional optics. 

- Use of binary optics with conventional optics to form an 
equivalent asphere or achromatized element. 

- Potential for low cost replication. 

- A reduction in the optical elements required thereby 
resulting in a weight savings of approximately 33%. 

- Coherent laser array beam addition to form a single 
powerful laser from many low power lasers. 

- Laser beam profile shaping for more uniform intensity 
distribution for applications like laser radar. 

- Laser multiplexer for use in laser radar. 

- IR broadband imaging systems and IR filters for use in 
optically pumping far IR lasers. 

- Solar cell concentrators to improve efficiency and reduce 
vulnerability to antisatellite threats. 

- Flexibility in optical system design to generate arbitrary 
compact optical shapes which produce the desired optical 
wave front. 

Other applications include helmet-mounted displays, 
uncooled thermal imagers, RLG readout optics, antiref lective 
avionic displays, optical storage disc readout, laser beam 
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steering, and advanced electro-optic image processing at the 
detector (i.e., amacronics) . 

B. NEAR-ANGLE SCATTERING ANALYSIS AND RECOMMENDATIONS 
1. Near-Angle Scattering Analysis 

Results are very promising for using a scalar theory 
to predict the near-angle scattering from binary optics. The 
advantages of using a scalar theory are: 

- The scalar theory is relatively straight-forward. Full 
EM theory is not required (i.e., polarization and phase 
angle of the light and boundary conditions are not 
required) to give an adequate approximation of scattering 
performance . 

- The scalar theory is not computer intensive and may be 
implemented on an IBM personal computer, resulting in 
significant cost savings in comparison to computer 
simulations requiring full EM theory. 

The shape of the theoretical curves closely follows 
the shape of the measured lobes in the central angular 
regions. The most disagreement with the scalar theory 
occurred at extremely small angles. Overall, experimental 
results agreed well with theoretical results. However, values 
for the autocorrelation length and surface roughness were much 
higher than the Talystep measurements indicated, which may be 
due to the short scan length limitation of the Talystep 
instrument. Longer autocorrelation lengths can be justified 
using a Wyco surface prof ilometer . 

Unfortunately, the binary optic used had known 
manufacturing defects (i.e., it was etched to an improper 
depth) and time constraints precluded obtaining a better 
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specimen. Keep in mind that this is a very new technology, 
hence an imperfect sample was generated in the process of 
"working out the bugs" with the manufacturing process. 
Consequently, this should not be considered detrimental to 
this promising technology nor as a typical occurrence. 

2 . Recommendations 

Since the scalar scattering theory developed depends 
on the autocorrelation length due to micro-roughness on the 
surface, future work is needed to determine how to best 
measure the autocorrelation length using either the Talystep 
or Wyco prof ilometer . Additional research could modify the 
theory by introducing a random offset to the periodicity of 
the surface structure (i.e., binary profile). 

DoD should implement binary optics in future electro- 
optic systems or upgrades to current military and space 
systems. Binary optics is a worthwhile area for research 
dollars. The ability to improve optical system performance 
with potential significant cost savings as well as perform 
functions impossible with conventional lenses, makes binary- 
optics technology a true revolution that will change the 
electro-optics industry. 
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APPENDIX A 



DEVELOPMENT OF SCALAR SCATTERING THEORY FOR BINARY OPTICS 



Equation A.l represents Fraunhofer diffraction from a 
binary optic and is evaluated in this appendix. 



U(x'.y') 



cxp(ik’/.) 

iXz 
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exp 



2k\ p • p ’ 
Xz 
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where 



Z exp 

p 



/ 2nipx \ 

\ D / 



exp 



/ 2ni?(|»(n- 1)\ 

l / 



= Fourier series expansion of phase changes due 
to the binary height variations of surface, 

= phase change due to the random roughness, 



£ (p) = random roughness on surface, 
n = index of refraction (constant) , 

D = period of binary optic (constant) , 

a p = Fourier coefficients (calculated for binary surface 
profile) , 

\ = wavelength (constant) , 

b = radius of the beam width (constant) , and 




= Gaussian laser beam with unit amplitude. 



For convenience in notation, define p = x x + y y, where 



x and y are unit vectors along the x and y axes respectively 
in the source plane, and p = x x' + y y', where x' and y' are 
unit vectors along the x 1 and y' axes respectively in the image 
plane . 

The point spread function (PSF) is equal to the square of 
Equation A.l: 
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Next, a change of variables is made. Let r = p - p" , 
t x = x - x", and r y = y - y". Making the change of variables 
in Equation A. 2 results in: 




Ensemble averaging both sides of Equation A. 3 yields, for 
Gaussian random variables: 1 




2k i (n - 1) 5(x+p")-$(p") 





(A. 4) 



where 

< > denotes the ensemble average, 

<5 = rms surface roughness, and 
G ( r ) = autocorrelation function. 

Substituting Equation A. 3 into the ensemble average of 
Equation A. 2 yields: 



^ung, C. C. and Holzer, J. A., "Scattering of Electromagnetic 
Waves from a Rough Surface," Applied Physics Letters, Vol. 28, No. 
8, pp. 429-431, 1976. 
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(A. 5) 



Treating the random roughness as a perturbation on the 
binary profile, first-order perturbation theory is applied for 
near angle scatter since <S << X. 
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Substituting A. 6 into A. 5 resullts in: 
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Xz 




dx-«iy"ZI»A 




-f oo 



/ 27tipx x \ 



1 - 



v q 



2rc(n - 1)5 



■G(i) 



j - 2ni p ’ • t 

-_ ; ox p 



/ r ' 2 \ 

TV - / (A - 7 > 



First, expand the dr dr integral in Equation A. 7. 
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(A. 8) 



The first integral in A. 8 is the specular component and 
the second integral is the scatter from the surface roughness. 
For the second integral in Equation A. 8, assume the form of 
the autocorrelation function is exponential: 
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(A. 9) 




where a = autocorrelation length. 

Substituting Equation A. 9 into the second integral of 
Equation A. 8 results in: 



w> 



roughness “ 





(A. 10) 



The following approximation can be made since the 
exp{ - 1 t | /a) falls off much more rapidly than the 

exp{-|p" + r| 2 /b 2 } for a < b: 




To show Equation A. 11 is a valid approximation, typical values 
are substituted into the left hand side of Equation A. 11 
(a = 100 /j m, b = 5000 \i m, and r = 2000 /xm) : 
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(A. 13) 



(A. 14) 
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Since exp{-|T| /a) is equal to values from 1 to 0 while 
exp{-|p" + r| 2 /b 2 } ~ exp{ - | p" | 2 /b 2 } , the approximation can be 
made that exp{-|p" + r| 2 /b 2 } ~ exp { - 1 p " 1 2 /b 2 ) for a < b. 

Substituting A. 11 into A. 7 using only terms associated with 
scattering due to the roughness results in: 




Evaluating the second integral in Equation A. 15 results 

in: 
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a = • 



P> = 



2rcm 

D 

2tix' 



Xz 



P 

Py Xz 



(A. 17) 
(A. 18) 
(A. 19) 



Evaluating the first part of the integral in A. 15, results 



m: 



b 2 re v v * / - b 2 tc 2 {p - q ) 2 

— — jl IwH 7 ~ 

2 (Xz) p ■« l 



2D 



(A. 20) 



Since exp( 



• b 2 7i 2 |p - q) 2 
2D 2 i 



0 for p f q, this simplifies to: 
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(A. 21) 



b k vi i 2 

ilw 

2(Xz) P 

Using the results from A. 16 and A. 21, the integral in A. 15 
simplifies to: 



<JU| 2 >; 



roughness — 



b 2 7i 



2(Xz) p 



2n(n - 1)5 



n 2 



2 no 





1 „ ,2 2 


2 


1 + 


[(« - P J + P yj 


a 



i 3 A 



(A. 22) 



A further simplification of A. 22 is possible if we examine 
the scatter in the plane (i.e., y* = 0 due to horizontal 
scan) . 



1 1 Xoughncss - y,| a p| 
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2 2 2 
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2n(n - 1)5 
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i 



1 + 
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(a-p x ) o 



3 / 2 



(A. 23) 



Next, the specular scatter from the first integral in A. 8 
is evaluated resulting in: 
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2n(n - 1)5 


2 


i - 


X 





b‘71 cxp(27ti.\"(^— exp/- b 2 K 2 l— - £) \cxp/2jii£^\ exo/ 



\ d i/^\ i r x T irr xp \ 



cxpJ 2ni¥—¥— J exp: - b 2 k 2 (^ 



,, 2 ) 



Xzj j 
(A. 24) 



Substituting A. 24 into A. 7 using only terms associated with 
the specular scatter results in: 



W> 



specular 



b 2 7i 






i - 



2n(n- 1)8 



// 



+ oo 



P q ' 1 \ 



• X 



ii 2 



exp' 



’ y ' exp/ 27tix" exp(-b 2 n 2 (^---^| \ exp ( 271 i exp(- b 2 n 2 (^-j 



b 



\ Uz D|/ 



Xz 



\ ^z / 



Xzj 



(A. 25) 
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Evaluating Equation A. 25 yields: 



W> 



specular 




2rc(n * 1)8 



/ r ‘ 

expj-h 71 
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(A. 26) 



If p = q 



<juf> 



specular = 



/.2 ^ 
bn 



Xz 
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X|a p | 2 expj-2b 2 7t 2 
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2 , , v 2 
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(A. 27) 



A further simplification of A. 27 is possible if we 
consider scatter only in the horizontal plane (i.e., y' = 0) . 



W> 



specular ” 



b 2 71 

Xz 



i - 



2n(n - 1)5 






(A. 28) 



The integral in Equation A. 7 has been evaluated and 
simplified using the results from A. 23 and A. 28. 
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(A. 29) 
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(A. 30) 



(A. 31) 



(A. 32) 
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APPENDIX B 



DETERMINATION OF FOURIER COEFFICIENTS FOR BINARY OPTICS 



Equation B.l represents the integrals used to determine 
the Fourier coefficients from the binary optic's surface 
profile (refer to Figure 5.10 in Chapter V). The integrals 
in B.l are evaluated in this appendix. 



J \ >• / 

" - on 




(B.l) 



where 

n = index of refraction, 

h(x) = surface profile as a function of x, and 
a p = Fourier coefficients. 

Using the orthogonality relation, Equation B.l becomes: 




exp 



j - 27tiqx ' 

i L - 



dx 



* +L /2 





a p exp 



/2m(p-q)x) 

\~L~j 



dx 



(B.2) 



where L is the period of the step structure on the binary 
optic . 

Evaluating B.2 results in: 




(B.3) 



if p = q, then 
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( B . 4 ) 



sinWjMl)) _ , 
njp^jl 
L 



if p t q, then 



sin(jt(|)-(])| _ 

7i(p-q) 

L ( B . 5 ) 



Consequently, B.3 simplifies to: 

^2 

CX p/ ^i|»W( n ^ Vcxp; 



/h 



/ 1 2nipx\ x = L; 



\ L / 



( B . 6 ) 



Next, the left side of the integral in B.6 is expanded: 



x x 
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exp 
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4* X 2 
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-X, 



cxp j^p idx 



+ exp/^Bll) 



• L /2 



+ X *■ 






(B.7) 



where H is the height of binary optic step. Evaluating B.7 
and simplifying yields: 



expJMSiUn 



. l- n P x 2 
sin 



(2rcpx,\ . |2npx,\ . (2npx 2 

- sin(np) - sinl — j- — + sinl — j- — I -sinl 



•JtP 

L 



( B . 8 ) 



Setting B.8 equal to La p and dividing through by L yields: 



a p = 






. /27tpx 2 \ . I2npx , 

sin[ — “ — - sin(7tp) - sinl — — — 



(2npx,} . [2jtpx 2 

+ sinl— j— J - 



Ttp 



( B . 9 ) 
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Collecting terms in B.9 and simplifying for p t 0 results in: 



a i> = 






/27ipxj\ . /27tpx 2 
s,n ‘ L ] ’ L~ 



-Tip 



Now consider the case for p = 0, Equation B.7 becomes: 






-X 



\ r + x ' 

dx + exp^ — - > dx 

X 2 \ X p-x, 



f 27 tiH(n-l) 



.+ X, 



•I dx + exp/ — - 

J+Xy \ X 



r 

J+x 



Vi 



dx 



Evaluating B.ll and simplifying: 



, , / 2niM(n-l)\ « lL 

a 0 L= 2 (x 2 - x j) + exp | t|2|y-x 2 +: 



B. 10) 



B. 11) 



B. 12) 
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APPENDIX C 



fortran program to model scalar scattering 

FROM A BINARY OPTIC 



l.T I -ARRY V. CIIIZKK 
NAVAL POSTGRADUATE SCHOOL 
29 MAY 1989 

FORTRAN PROGRAM TO CALCULATE TRANSMISSIVE NEAR ANGLE SCATTER FROM 
BINARY OPTIC USING SCALAR THEORY (FRESNEL DIFFRACTION INTEGRAL) 
FOR NORMAL INCIDENCE. 

PROGRAM ASSUMES GAUSSIAN FOR SPECULAR AND 
SCATTERED COMPONENT WITH EXPONENTIAL AUTOCORRELATION 

DECLARE VARIABLES 

INTEGER I . J , M , NPTS , T 



REAL X IN , ZM , PI , X 1 , X2 , II , N , L, LAMBDA , DELTA , Cl , C2 . ALPHA , BETA 

REAL IMPINT,A(51) . RL(51) , IM(51) , U2 (42) ,THETA(42) ,£V(A2) ,EV25(42) 

REAL SUM1 , SUM2 . INCR , CL, STEP . BU , CAUSEXP , C3 , C4 , PKD , PKE , PK , DIS P 

COMMON A ,RI. , IM , PI , II , N , I.,X1 , X2 , IAMISDA , T 



OPEN (1, FILE- 


’C : EXP1L5.DAT 1 ) 






OPEN (2, FILE- 


• D: THEORY. DAT' ) 






OPEN (5.FILE- 


■C: EXP1L25.DAT' ) 






DISP 


DISPLACEMENT CONSTANT FOR LOBES 






PKD 


PEAK OF CALCINATE!) DATA 






PKE 


PEAK OF EXPERIMENTAL DATA FOR 500 PINHOLE 




PK 


PEAK OF EXPERIMENTAL DATA FOR 25 PINHOLE 




EV 


EXPERIMENTAL VALUES FROM 500 MICRON PINHOLE 




EV25 


EXPERIMENTAL VALUES FROM 25 MICRON PINHOLE 




SUMl 


CONSTANT 






SUM2 


CONSTANT 






GAUSEXP 


EXPONENT FOR GAUSSIAN FUNCTION 






BW 


BEAM WIDTH AT OPTIC 






NPTS 


NUMBER OF DATA PTS TO BE GENERATED 






INCR 


INCREMENT FOR DATA PTS IN INCHES 






A(4) 


FOURIER COEFFICIENTS 






RL(4) 


REAL PART OF FOURIER COEFF 






IM(4 ) 


IMAGINARY PART OF FOURIER COEFF 






L 


PERIOD OF BINARY OPTIC 






XI 


TRANSITION PT FOR HEIGHT ON SURFACE OF 


OPTIC 


IN METERS 


X2 


TRANSITION PT FOR HEIGHT ON SURFACE OF 


OPTIC 


IN METERS 


H 


HEIGHT OF BINARY OPTIC IN METERS 






ZM 


DISTANCE FROM BINARY OPTIC TO DETECTOR 


IN METERS 


ZIN 


DISTANCE FROM BINARY OPTIC TO DETECTOR 


IN INCHES 


LAMBDA 


WAVELENGTH OF HENE LASER 






N 


INDEX OF REFRACTION FOR BINARY OPTIC 






DELTA 


RMS SURFACE ROUGHNESS OF OPTIC IN METERS 




U2 


INTENSITY OF SCATTERED LIGHT 






Cl 


CONSTANT 






C2 


CONSTANT 






C3 


CONSTANT 






C4 


CONSTANT 






ALPHA 


CONSTANT 






BETA 


CONSTANT 






THETA 


SCATTER ANGLE IN RADIANS 
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C 

C 

C 



I 

J 

M 



PI 



STEP 



COUNTER 

COUNTER 

COUNTER 



PI RADIANS 

INCREMENT FOR ANGULAR RESOLUTION IN INCHES 



INITIALIZE VARIABLES 

DISP - 0 . 3164E- 2 
H - 6554E- 10 
L - 200E- 6 
XI - 42.92E-6 
X2 - 86.5E-6 
N - 1.46 

LAMBDA - 6328E- 10 
DELTA - 80.0E-10 
PI - 3.1416 
CL - 175.0E-6 
ZIN - 71.6 
ZM - 1.819 
INCR - 0.005 
NPTS - 42 
BW - IE- 2 
PKD - 0 
PKE - 0 
PK - 0 

INITIALIZE ARRAYS 

DO 100 I - 1 , NPTS 

U2(I) = 0 
EV(I ) - 0 
EV25 (I) = 0 
THETA(I) ~ 0 

100 CONTINUE 

DO 200 J - 1,51 



200 CONTINUE 

CALCULATE ANGLE IN RADIANS 
DO 300 I - 1 , NPTS 



STEP ~ I * INCR 

THETA(I) - ALOGIO (STEP/ZIN) 



A( J ) - 0 
RL(J ) - 0 
IM(J ) - 0 



300 CONTINUE 



C 



CALCULATE CONSTANTS 
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Cl - (1 - ((2.0*PI*(N-1 )*DELTA) /LAMBDA )**2 )* ( l/( LAMBDA*ZM) )**2 
C2 - (((2.0*PI*(N-l)*DELTA)/LAMBDA)**2)*((BW/2)/(LAMBDA*ZM))**2 
C3 - ( (BW/2)**4)*(PI**2) 

C4 - Cl * C3 



C 



C 

c 

c 

c 

c 

c 

c 

c 

c 

c 



/ 

c 

c 

700 

c 

c 

c 

c 

c 



c 

c 

c 

c 

c 

c 

c 

c 

c 



CALCULATE INTENSITIES FOR SCATTERING ANGLES 
DO 600 I - 0 , NPTS- 1 

DO 700 M - -25,25,1 
T - 26 + M 

STEP - ( ( I*INCR)/ZIN) + (1 . 0*DISP) 

CALL FRCF(M) 

GAUSEXP - -2.0*(((PI*BW)/2)**2)*((STEP/LAMBDA 
- (M/L)))**2 

SUM1 - A(T)*C4*EXP(GAUSEXP) 

ALPHA - ( 2 . 0*PI*M)/L 

BETA - (2.0*PI*STEP)/LAMBDA 

SUM2 - (A(T)*( (PI**2)*(CL**2)*C2) )/ 

(1 + (CL* (ALPHA - BETA) )**2)**1 . 5 

U2 (1+1) - SUMl + SUM 2 + U2(I+1) 

CONTINUE 

READ (5 , *) EV25 ( 1+1 ) 

READ (1 , *) EV(I+1) 

NORMALIZE LOG OF DATA 
IF (I .NE. 0) THEN 
PKD - U2 ( 1 ) 

TAKE LOG AND NORMALIZE EXP DATA 
PKE - EV25 (1)*(244 . 0) 

PK - EV25 (1 ) 

U2(I+1) - U2(I+1)/PKD 
U2(I+1) - ALOGIO (U2(I+1)) 

EV(I+1) - EV(I+1)/PKE 
EV(I+1) - ALOGIO (EV ( 1+1 ) ) 
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C 

C 



/«/* 



c 

c 



EV25(I+1) - EV25 ( 1+ 1)/PK 
EV25(I+1) - A LOG 10 ( EV25 ( 1+1 ) ) 

WRITE(2 , 4A ) THETA (1+1) , EV25(I+1) , EV(I+1) ,U2(I+1) 
FORMAT ( 1 X , ' ()• , E 1 0 . 3 , 1 X . K1 0 . 3 , 1 X , K1 0 . 3 . 1 X , FI 0 . 1 ) 



ENDIF 



600 CONTINUE 



STOP 
END 

SUBROUTINE FRCF CALCULATES MAGNITUDE SQUARED OF FOURIER COEFFICIENTS 
FOR SURFACE OF BINARY OPTIC 

SUBROUTINE FRCF(M) 

REAL K1 f K2,K3,RL(51) ,IM(51) ,A(51) , LAMBDA , H , N , PI , L, XI , X2 
INTEGER M,T 

COMMON A , RL , IM , PI , H , N , L , XI , X2 , LAMBDA , T 

K1 CONSTANT 

K2 CONSTANT 

K3 CONSTANT 

CALCULATE CONSTANTS 

K1 - ( 2 . 0*PI*H* (N- 1) ) /LAMBDA 
K2 - (2 . 0/L)*( (L/2) - X2 + XI) 

K3 = SIN( ( 2 . 0*PI*M*X1 ) /L) - SIN( (2 . 0*PI*M*X2)/L) 

C 

IF (M .EQ. 0) THEN 
C 

RL(T) - (2.0*(X2-X1)/L) + (K2 * COS(Kl)) 

IM(T) - K2 * SIN(Kl) 

A(T) - (RL(T)**2 + IM(T)**2) 



RL(T) - ((1.0 - COS(Kl))*K3)/(-PI*M) 
IM (T) - (SIN(K1)*K3)/(PI*M) 

A(T) - (RL(T)**2 + IM(T)**2) 

C 

ENDIF 

C 

RETURN 

C 

END 
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